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CHAPTER  1 


INTRODUCTION 


1.1.  Problem  Introduced 

There  has  been  increasing  interest  in  recent  times  in  the  study 
of  unsteady  viscous  flows.  This  is  not  surprising  since,  to  a 
certain  degree,  unsteady  viscous  effects  are  present  in  the  majority 
of  practical  flow  situations.  Viscous  effects  play  either  the 
leading  or  a  secondary  role  in  such  varied  fields  as  biofluid  flows, 
turbomachinery  flows,  aircraft  flutter,  helicopter  rotor-blade  flows, 
etc. 

The  particular  case  of  a  boundary  layer  subjected  to  a  free- 
stream  velocity  which  is  forced  to  change  sinusoidally  with  time 
around  a  nonzero  mean  is  of  primary  theoretical  and  practical 
importance.  Its  study  received  pioneering  contributions  from 
Lighthill  (1954)  and  Karlsson  (1959).  The  former  introduced  an 
analysis  and  the  first  (asymptotic)  solutions  for  laminar  flow,  which 
set  the  basis  for  future  theoretical  investigations  on  the  subject. 
Karlsson,  in  turn,  performed  the  first  series  of  experiments  in  the 
turbulent  regime,  covering  a  wide  range  of  cases  (varying  the 
frequency  and  amplitude  of  the  imposed  free-stream  oscillation)  which 
have  not  been  matched  by  any  other  study,  either  in  number  or  in  the 
range  of  experimental  conditions. 


Since  then,  broadly  three  different  lines  of  approach  have  been 
taken  to  study  the  phenomenon  of  oscillatory  boundary  layers.  One  is 
the  theoretical  analysis,  which  has  been  mainly  confined  to  laminar- 
flow  situations.  In  the  restricted  cases  it  can  deal  with,  it 

provides  an  understanding  of  the  interplay  among  the  different 

physical  mechanisms  and  the  resulting  effects.  The  works  of 
Ackerberg  and  Phillips  (1972)  and  Pedley  (1972)  belong  to  this 
category.  In  turbulent  flow,  however,  the  theoretical  analysis  is 
difficult  because  of  the  well-known  closure  problem. 

Another  line  of  approach  is  the  numerical  solution  of  the 

differential  equations.  In  laminar  flow,  this  allows  one  to  extend 
the  analytical  solutions  beyond  their  range  of  validity,  sometimes 
bridging  the  gaps  between  asymptotic  solutions.  Reference  can  be 

made  to  the  works  of  Tsahalis  and  Telionis  (1974)  and  McCroskey  and 
Phi  11 i pe  (1975)  in  this  connection.  When  dealing  with  turbulent 
flows,  the  use  of  a  turbulence  model  is  necessary  in  order  to  close 

the  system  of  equations.  This  limits  considerably  the  generality  of 

the  solution.  In  addition,  in  unsteady  flow,  the  turbulence  models 
used  so  far  are  unmodified  versions  of  those  developed  for  steady 
flow.  Therefore,  the  predictions  are  expected  to  fail  for  relatively 

large  frequencies  of  oscillation.  Work  in  this  area  has  been  done 

by,  among  others,  Nash  and  Patel  (1975),  Cousteix  et  al.  (1977), 
Orlandi  (1981)  and  Menendez  and  Ramaprian  (1982,  1983a,  1983b). 
Integral  methods  have  also  been  developed  to  obtain  approximate 
solutions  [Cousteix  et  al.  (1977),  Lyrio  and  Ferziger  (1983)]. 


The  third  line  of  approach  to  the  study  of  unsteady  boundary 
layers  is  through  experiments.  With  regard  to  laminar  flows  they 
have  served  primarily  to  corroborate  the  theoretical  results.  Hill 
and  Stenning  (1960)  have  provided  some  of  the  few  known  experimental 
data.  In  the  case  of  turbulent  flows,  on  the  contrary,  much  of  the 
current  understanding  of  the  phenomenon  comes  from  the  experimental 
studies.  These  include  the  works  of  Patel  (1977),  Schachenmann  and 
Rockwell  (1976),  Cousteix  et  al.  (1977,  1981),  Jarayaman  et  al. 
(1982),  and  Simpson  et  al.  (1983a,  1983b).  Related  investigations  on 
duct  (pipe  or  channel)  flows  have  been  performed  by  Binder  and  Kueny 
(1981),  Kobashi  and  Hayakawa  (1981)  and  Tu  and  Ramaprian  (1983a, 
1983b).  In  a  recent  book,  Telionis  (1981)  summarizes  much  of  the 
current  knowledge  about  unsteady  boundary  layers. 

1.2.  State  of  the  Art  in  Unsteady 
Turbulent  Boundary  Layers 

From  the  available  experimental  information,  the  following 
picture  of  the  phenomenon  of  oscillatory  turbulent  boundary  layers 
can  be  constructed: 

(i )  The  time-mean  velocity  distribution  is  practically 

coincident  with  that  in  steady-flow  at  the  time-mean  free-stream 
velocity  [Karlsson  (1959),  Cousteix  et  al.  (1977),  Jarayaman  et  al. 
(1982)].  Small  deviations  from  this  behavior  were  observed  by 

Ramaprian  and  Tu  (1982)  in  pipe  flow. 


(ii)  The  oscillatory  velocity  component  depends  strongly  on  the 
forcing  frequency.  For  low  frequencies,  the  unsteady  effects  spread 
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over  the  outer  part  of  the  boundary  layer.  For  larger  frequencies, 
the  outer  part  is  subjected  to  rigid-body  oscillations,  and  the 
unsteady  viscous  effects  are  concentrated  in  regions  increasingly 
closer  to  the  wall  [Karlsson  (1959),  Jarayaman  et  al.  (1982)]. 
However,  unsteady  viscous  effects  in  turbulent  flow  extend  over  a 
much  wider  region  than  in  a  hypothetical  laminar  flow  with  the  same 
boundary-layer  thickness  and  at  the  same  frequency  of  oscillation. 
Using  eddy-vi scosi ty  arguments,  Ramaprian  and  Tu  (1982)  identified 
the  frequency  parameter  036/u^  as  the  one  characterizing  the  extent  of 
unsteady  viscous  effects  in  the  outer  layer.  They  also  used  this 
parameter  to  define  five  different  frequency  regimes.  However,  their 
analysis  is  not  rigorous  but  heuristic.  Furthermore,  it  does  not 
properly  account  for  the  viscous  effects  very  near  the  wall  and  hence 
is  unable  to  characterize  the  flow  completely,  especially  at  high 
frequencies. 

(iii)  The  amplitude  of  the  oscillatory  velocity  shows  an 
overshoot,  which  is  especially  large  for  adverse  time-mean  pressure- 
gradient  flows  at  low  frequencies  [Cousteix  et  al.  (1981),  Jarayaman 
et  al.  (1982)].  For  high  frequencies,  this  overshoot  is  small 
irrespective  of  the  free-stream  conditions.  Close  to  the  wall  and, 
at  least,  for  relatively  large  frequencies,  the  oscillatory  velocity 
leads  the  free-stream  velocity  [Karlsson  (1959)].  At  high 
frequencies,  the  oscillatory  motion  seems  to  approach  the  laminar 
(Stokes)  solution  given  by  Lighthill  (1954)  [Binder  and  Kueny  (1981), 
Jarayaman  et  al.  (1982)]. 
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(iv)  The  turbulence  properties  are  modulated  with  the  forcing 
frequency.  Some  of  the  existing  data  suggest  that  their  time-mean 
distributions  coincide  with  the  steady-flow  distributions  at  the 
time-mean  free-stream  velocity  [Cousteix  et  al.  (1977),  Jarayainan  et 
al.  (1982)].  Others  [Ramaprian  and  Tu  (1982)]  do  not  seem  to  suggest 
this  trend. 

(v)  The  oscillatory  component  of  the  turbulence  properties  has 
a  strong  dependence  on  the  forcing  frequency.  For  low  frequencies, 
significant  amplitude  and  phase  variations  occur  all  across  the 
boundary  layer.  For  higher  frequencies,  these  variations  are 
confined  to  the  near-wall  region,  and  the  turbulence  oscillations  are 
practically  supressed  in  the  outer  layer  [Jarayaman  et  al.  (1982)]. 

(vi )  There  is  controversy  about  the  interaction  between  the 
forced  oscillation  and  the  turbulent  fluctuations.  It  still  remains 
to  be  answered  if  there  is  any  interaction;  and,  if  so,  at  what 
frequency  and  what  its  nature  is. 

(vii)  No  direct  wal  1 -shear-stress  measurements  have  been 
performed.  Cousteix  et  al.  (1981)  utilized  the  conventional 
procedure  developed  for  steady  flow  to  obtain  the  wall  shear  stress 
distribution  in  a  cycle.  This  consists  of  assuming  the  validity  of 
the  "logarithmic  law"  for  the  ensemble-averaged  velocity  profile  and 
using  the  "Clauser  plots".  However,  it  is  not  clear  whether  this 
procedure  is  correct.  Moreover,  for  high  frequencies  the  logarithmic 
region  may  disappear  for  part  of  the  cycle  [Jarayaman  et  al. 
(1982)].  Ramaprian  and  Tu  (1982)  performed  direct  wal 1 -shear-stress 
measurements  in  pipe  flow.  They  found  that  the  time-mean  value  was 
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slightly  larger  than  the  steady  value  at  the  time-mean  free-stream 
velocity. 

In  an  admirable  effort,  Carr  (1981b)  is  compiling  all  the 
available  experimental  data  on  unsteady  turbulent  boundary  layers. 
This  has  provided  him  with  the  opportunity  of  establishing  some 
correlations  among  the  different  experimental  conditions,  which  may 
bring  more  order  into  the  ongoing  research.  Figures  1.1  and  1.2, 
taken  from  Carr  (1981a),  show  two  different  ways  of  correlating 
experiments  with  zero  (ZPG)  and  adverse  (APG)  time-mean  pressure- 
gradient  boundary  layers  and  with  pipe  flows.  Both  are  in  terms  of 
the  frequency  and  amplitude  of  oscillation  and  the  extent  of  unsteady 
viscous  effects.  In  Fig.  1.1,  the  frequency  is  presented  in  terms  of 
a  Strouhal  number  based  on  the  local  time-mean  boundary  layer 
thickness  and  free-stream  velocity.  In  Fig.  1.2  f  is  normalized  by  a 
critical  value  associated  with  the  local  turbulent  burst  frequency. 
These  figures  suggest  that  there  is  still  a  need  for  experiments  on 
boundary  layers  at  zero  time-mean  pressure  gradient,  especially  at 
high  frequencies  and  large  amplitudes  of  oscillation. 


1.3.  Objectives  of  the  Present  Research 


The  objectives  of  the  present  work  were  twofold.  First, 
detailed  and  high  quality  experimental  information  was  sought  for  a 
turbulent  boundary  layer  in  a  zero  time-mean  pressure  gradient  with  a 
relatively  large  amplitude  of  oscillation.  Two  frequencies  were 
studied;  for  the  lower  one,  the  unsteady  viscous  effects  extend  over 


the  outer  la-'er,  while  they  are  confined  in  the  near-wall  region  for 


the  higher  frequency.  The  higher  frequency  is  comparable  to  the 
turbulent  bursting  frequency  in  the  boundary  layer.  One  can, 
therefore,  expect  the  imposed  oscillation  to  interact  with  the 
turbulent  motions.  The  experiments  are  represented  in  Figs.  1.1  and 
1.2  by  means  of  lines,  which  cover  the  conditions  from  the  first  to 
the  last  measuring  station  along  the  tunnel  employed.  The  main 
differences  with  the  similar  experiments  performed  by  Karlsson  (1959) 
and  Cousteix  et  al.  (1977)  are: 

(i )  Water  was  used  instead  of  air  as  the  working  fluid.  Hence, 
the  amplitude  modulation  imposed  by  acoustic  effects  was  avoided  and 
a  long  test  section  could  be  used.  This  allowed  study  of  spatial 
history  effects  as  the  boundary  layer  evolves  through  transition  and 
develops  into  a  fully  turbulent  boundary  layer. 

(ii)  The  high  frequency  utilized,  when  properly  normalized,  was 
much  larger  than  those  used  by  Cousteix  et  al.  (see  Figs.  1.1  anti 
1.2).  Though  some  of  Karlsson's  frequencies  are  comparable  to  the 
present  ones  (or  even  much  larger  than  the  present  ones,  as  in  an 
experiment  not  represented  in  the  foregoing  figures),  he  was  not  able 
to  isolate  properly  the  turbulence  properties  due  to  limitations  in 
his  data  processing  system. 

(iii)  Laser  Doppler  Anemometry,  capable  of  nonintrusive 
measurements,  was  used  in  the  present  studies.  Also,  wall -shear- 
stress  measurements  were  made  using  a  surface-mounted  probe  rather 
than  inferring  it  indirectly  from  a  'law  of  the  wall'  of 
unestablished  validity. 


Some  similarities  between  the  present  experiments  and  those  by 
Jarayaman  et  al.  (1982)  are  evident.  However,  the  latter  correspond 
to  an  adverse  time-mean  pressure-gradient  flow.  In  addition,  in  the 
present  experiments,  two  velocity  components  were  measured  instead  of 
one. 

The  second  objective  of  the  present  research  was  to  provide, 
following  an  analytical  approach,  a  theoretical  framework  to  classify 
and  characterize  different  unsteady  turbulent-flow  regimes  for 
nonseparating  boundary  layers  in  an  arbitrary  (time-mean)  pressure 
gradient,  as  well  as  for  fully  developed  pipe/channel  flows.  A 
singular  perturbation  technique  was  used,  and  led  to  asymptotic 
similarity  laws  for  the  different  frequency  regimes.  These  laws  can 
be  considered  as  extensions  to  unsteady  flows  of  the  classical 
steady-flow  laws,  such  as  the  "law-of-the-wall "  for  the  inner  layer, 
the  "velocity -defect  law"  for  the  outer  layer,  and  the  "logarithmic 
law"  for  the  overlap  region. 


1.4.  Layout  of  the  Report 


The  experimental  apparatus  and  procedures  are  described  in 


Chapter  2.  The  results  of  the  experiments  are  presented  and 
discussed  in  Chapter  3.  The  asymptotic  analysis  of  unsteady 


turbulent  boundary  layers  and  pipe/channel  flows  is  developed  in 


Chapter  4.  Chapter  5  presents  an  examination  of  the  present  as  well 


as  earlier  experimental  data  in  the  light  of  this  analysis.  The 
major  conclusions  of  the  present  work  are  summarized  in  Chapter  6. 
The  detailed  experimental  data  are  too  voluminous  to  be  presented  in 
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Figure  1.2.  Correlation  between  amplitude  and  relative 

frequency  [Carr  (1981a)].  - ,  present  experi¬ 

ments  . 


CHAPTER  2 


EXPERIMENTAL  DETAILS 

2.1.  General 

The  object  of  the  experimental  investigation  was  to  generate  and 
study  a  zero  time-mean  pressure  gradient  boundary-! ayer  flow,  with  a 
free-stream  velocity  given  by 

<U  >  =  U  ( 1+e  sin  u>t)  (2.1) 

where  U0,  e  and  a>  =  2irf  are  constant.  A  large  water  tunnel  was  built 
for  this  purpose.  The  use  of  water  as  the  working  fluid  has  the 
following  advantages: 

(i)  It  is  possible  to  study  oscillation  frequencies  that 

approach  the  turbulent  frequencies. 

(ii)  There  is  a  significant  separation  between  the 

oscillation  and  turbulent  frequencies  on  one  hand,  and 

the  acoustic  frequencies  on  the  other. 

(iii)  No  seeding  is  necessary  for  use  of  Laser  Doppler 


Anemometer 
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In  addition,  there  exists  a  good  deal  of  experience  at  1 1  HR  in 
building  such  facilities  for  steady  flow. 

The  main  instrumentation  consisted  of  a  Laser  Doppler  Anemometer 
(LDA)  to  measure  the  velocity,  and  a  heat-flux  gage  (HFG)  operated  by 
constant-temperature  hot-wire  anemometer  to  determine  the  skin 
friction.  The  data  were  digitally  acquired  and  processed  on  an 
HP/1000  minicomputer.  More  details  are  presented  below. 

2.2.  Apparatus 

A  unique  water  tunnel  was  built  especially  for  this  study.  Fig. 
2.1  shows  a  layout  of  the  apparatus.  It  works  under  a  constant  head 
of  about  8  m.  The  test  section  is  2.4  m  long,  with  a  rectangular 
cross-section  50  cm  in  width  and  22.5  cin  in  height.  Its  bottom,  used 
as  the  test  surface  for  the  boundary  layer  studies,  is  a  smooth  brass 
plate.  Transition  was  promoted  by  means  of  a  strip  of  14-mesh  brass 
screen  (15  cm  in  width)  glued  to  the  test  surface  at  a  distance  (x) 
of  about  22.5  cm  from  the  end  of  the  contraction.  The  side  walls  of 
the  tunnel  are  9.5  mm  thick  steel  plates,  with  five  windows  located 
at,  approximately,  x  =  48,  69,  90,  142  and  203  cm.  The  windows  are 
made  of  Plexiglas.  Originally,  the  side  walls  were  tempered  glass, 
but  it  was  found  that  they  depolarized  the  laser  beams,  so  they  had 
to  be  removed.  The  top  wall  of  the  test  section  is  a  Plexiglas  sheet 
37.5  mm  thick.  In  addition,  a  6.4  mm  thick  flexible  Plexiglas  sheet 
was  used  as  a  'false'  wall.  By  adjusting  its  position  to  correct  for 
the  displacement  thickness  of  the  boundary  layers  along  the  walls  a 
nearly  zero  (time-mean)  longitudinal  pressure  gradient  was  obtained 
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in  the  test  section.  For  the  0.5-Hz  experiments  (see  below),  the 
false  wall  was  reinforced  with  transverse  strips  6.4  to  9.5  mm  thick, 
spaced  approximately  23  cm  apart.  Other  details  of  the  tunnel  can  be 
seen  in  Fig.  4.1.  Note,  in  particular,  that  no  screens  were  used  in 
the  tunnel.  This  was  to  avoid  clogging  problems,  likely  to  arise 
from  the  use  of  unfiltered  recirculated  water.  The  tunnel  ends  in  a 
steel  cylinder  (40  cm  diameter  x  60  cm  long)  whose  downstream  end  is 
closed.  The  cylinder  has  two  longitudinal  rectangular  slots,  61  x 
2.54  cm,  located  180  degrees  apart,  from  which  the  water  exits.  The 
area  of  opening  of  the  slots  is  varied  by  a  rotating  profiled  sleeve 
driven  by  a  3  H.P.  geared  D.C.  motor  whose  speed  can  be  regulated  to 
within  1/4  percent.  Each  complete  rotation  of  the  sleeve  corresponds 
to  two  oscillation  cycles.  The  sleeve  profile  was  contoured  to 
produce  a  sinusoidally  varying  free-stream  velocity  in  the  test 
section  at  the  desired  frequency.  The  details  of  the  design  are 
given  by  Ramaprian  and  Tu  (1982).  Note  that  a  different  profile 
should  be  used,  in  principle,  for  each  desired  frequency  (and 
amplitude)  of  oscillation.  However,  it  was  found  that,  with  a  given 
profile,  the  frequency  can  be  changed  appreciably  from  the  design 
value  without  producing  a  large  change  in  the  amplitude  or  any 
significant  distortion  in  the  velocity  wave  form. 

2.3.  Instrumentation 

Two  types  of  instruments  were  used;  a  two-component  LDA,  and  a 
flush-mounted  HFG.  The  first  measures  the  longitudinal  and  vertical 
velocity  components.  The  HFG  measures  the  shear  stress  at  the 


wall.  Each  of  these  is  described  in  some  detail  in  the  following 
subsections. 

2.3.1.  Laser  Doppler  Anemometer 

The  principle  of  operation  of  the  LDA  is  well  known.  Two  laser 
beams  intersect  each  other  producing  an  interference  fringe  pattern 
in  the  small  intersecting  region  ("measuring  volume").  When 
particles  suspended  in  the  fluid  cross  the  measuring  volume,  they 
scatter  light  with  a  frequency  proportional  to  the  particle  velocity, 
assumed  to  be  the  same  as  the  fluid  velocity.  The  detector  of  the 
f requency-shi ft  provides  a  signal  linearly  related  to  the  velocity. 
For  more  details  see,  for  example.  Drain  (1980). 

In  the  present  experiments,  a  two-component  LDA  was  used.  This 
system  works  with  three  laser  beams  which,  when  properly  polarized, 
produce  two  sets  of  interference  fringes,  allowing  the  simultaneous 
measurement  of  two  velocity  components.  The  system  consists  of  a 
SPECTRAPHYSICS  5  mW  He-Ne  Laser,  TSI  (Thermo  System  Inc.)  optics,  an 
acousto-optic  Bragg  cell  for  frequency  shifting  (which  allows 

directionality  distinction),  a  pair  of  photodetectors  (one  for  each 

signal)  and  a  pair  of  TSI  Model  1090  frequency  trackers,  which 
process  the  photodetector  signals.  The  output-voltage  signal  on  each 
channel  is  then  low-pass  filtered  by  an  active  filter  (Rockland, 
Model  852)  with  a  sharp  cut  off  (36  db/octave),  and  sampled 
simultaneously  by  the  computer.  A  sketch  of  the  optical  arrangement 
is  shown  in  Fig.  2.2.  For  more  details,  see  Ramaprian  and 

Chandrasekhara  (1983).  The  LDA  was  operated  in  the  forward  scatter 


inode.  A  beam  expander  was  used  to  reduce  the  size  of  the  focal 
volume  to  1.8  mm  x  0.1  mm.  The  transmitting  and  receiving  optics 
were  both  mounted  on  a  traverse,  which  could  be  moved  in  three 
spatial  directions. 

The  calibration  of  the  LDA  depends  on  the  frequency  of  the  laser 
beams  (fixed)  and  their  geometrical  configuration.  Once  the  latter 
is  set,  the  calibration  constant  is  fixed,  known,  and  independent  of 
the  experimental  conditions.  To  determine  the  distance  of  the 
measuring  volume  from  the  wall,  it  is  necessary  to  establish  the 
position  of  the  wall  ("zero1').  This  is  done  with  an  estimated  error 
of  the  order  of  the  laser  beam  thickness  at  the  focal  plane  of  the 
optical  system,  namely  0.1  mm. 


2.3.2.  Heat-Flux  Gage 

The  use  of  the  flush-mounted  HFG  in  steady  laminar  or  turbulent 
flow  is  based  on  the  analogy  between  heat  transfer  from  the  heated 
element  to  the  fluid  and  the  local  wall  shear  stress.  The  proposed 
relation  is  the  following: 


,1/3  =  AE3  +  B 

Mil 


(2.2) 


where  the  voltage  E,  read  from  the  anemometer,  is  related  to  the  wall 
heat  transfer,  and  the  constants  A  and  B  are  obtained  by  calibration 
[see  Bellhouse  and  Schultz  (1966)].  As  part  of  the  present  work, 
considerable  effort  was  directed  to  determine  the  extent  to  which 
this  technique  is  applicable  to  unsteady  flow  situations.  As  a 
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result,  a  more  general  relationship  between  the  two 
quantities,  i.  and  E,  has  been  developed,  namely  [Menendez  and 
Ramaprian  (1983c)] 


tw  =  (AE2  +  B)3  + 


C1  ^  +  c  A  dE^ 

(AE2+B)  dt  2  dt 


(2.3) 


where  the  constants  c^  and  C2  depend  on  the  fluid  properties  and  the 
effective  dimensions  of  the  heated  element.  Fig.  2.3,  taken  from  the 
above  reference,  shows  typically  the  improvement  in  results  for  phase 
and  amplitude  of  <tw>  obtained  by  the  use  of  Eq.  (2.3)  instead  of  Eq. 

(2.2) ,  under  different  operating  conditions.  The  results  are 
compared  in  each  case  with  theoretical  results  obtained  from  a 
numerical  solution  of  the  unsteacty  boundary  layer  equations.  Eq. 

(2.3)  reduces  practically  to  Eq.  (2.2)  when  the  frequency  of 
oscillation  is  small  enough  to  satisfy  the  conditions 


(2.4) 

(2.5) 


[For  details,  see  Menendez  and  Ramaprian  (1983c)].  These  conditions 

were  satisfied  in  the  present  experiments.  Hence,  Eq.  (2.2)  was  used 

for  the  determination  of  t, . 

w 

A  TSI  quartz-coated  film  probe  was  used,  and  operated  by  a  DISA 
constant-temperature  hot-wire  anemometer.  A  sketch  of  the  probe, 
mounted  on  the  wall,  is  shown  in  Fig.  2.4.  An  overheat  ratio  of 
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about  1.1  was  used.  Measurements  were  performed  at  four  stations;  x 
=  48  cm,  x  =  80  cm,  x  =  142  cm  and  x  =  203  cm.  The  signal  was  low- 
pass  filtered  before  being  sampled  by  the  computer.  Simultaneously, 
an  LDA  signal  corresponding  to  the  longitudinal  velocity  was  also 
sampled  as  a  reference.  The  probe  was  mounted  just  before  each 
experimental  run,  and  removed  and  cleaned  after  the  run.  The 
calibration  was  performed  in  steady  turbulent  flow,  for  which  the 
wall  shear  stress  is  known  from  a  Coles'  fit  to  the  velocity  profile 
(see  section  2.6.3).  The  calibration  depends  critically  on  both  the 
film  and  water  temperatures,  and  thus  could  change  considerably 
during  an  experiment.  Hence,  two  calibration  curves  were  obtained 
for  each  experiment,  one  immediately  before  and  one  immediately  after 
the  unsteady  run,  and  an  average  of  the  two  was  taken  as  the  final 
curve.  For  more  experimental  details,  reference  may  be  made  to 
Ramaprian  and  Tu  (1982). 

2.4.  Acquisition  and  Processing  of  the  Information 

An  HP/1000  data  acquisition  system  with  a  PRESTON  Analog-to 
Digital  Converter  was  used  for  both  data  acquisition  and 
processing.  A  triggering  mechanism  synchronized  with  the  sleeve 
rotation  provides  the  sampling  pulses  at  fixed  phase  positions.  This 
is  described  in  detail  in  Ramaprian  and  Tu  (1982).  An  improvement 
over  the  previous  way  of  operation  is  that  the  starting,  stopping, 
and  eventual  repetition  (in  case  of  missed  pulses)  of  the  sampling 
was  done  automati cal ly  by  the  computer.  The  triggering  mechanism 
gives  200  pulses  for  each  complete  revolution  of  the  sleeve,  which 
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corresponds  to  100  pulses  per  oscillation  cycle.  For  the  steady-flow 
experiments,  the  pulses  were  provided  by  a  square-wave  generator. 
For  each  pulse,  the  computer  sampled  simultaneously  the  two  signals 
(either  two  LDA  signals  or  one  HFG  and  one  LDA  signal)  received  from 
the  two  channels.  This  information  was  temporarily  stored  in  disk 
file  and  later  processed.  Both  the  raw  and  processed  data  were 
permanently  stored  on  digital  tape.  The  latter  type  of  files  are  in 
the  format  appropriate  for  the  AGARD  Data  Bank  of  L.W.  Carr  (1981b). 

The  software  for  data  acquisition  and  processing  of  single 
periodic  unsteady  signals,  previously  developed  at  1 1  HR ,  has  been 
extended  to  deal  with  two  simultaneous  signals.  It  has  also  been 
enhanced  with  new  features  which  allow  the  automatic  recognition  of 
spurious  signals  and  the  consequent  discarding  of  the  corresponding 
cycle  during  the  processing.  Spurious  signals  may  appear  due  to  LDA 
drop-outs  [see  Drain  (1980)]  or  the  temporary  locking  of  the  tracker 
to  a  secondary  signal  (for  example,  the  signal  corresponding  to  the 
other  channel). 

The  signals  are  processed  according  to  the  concept  of  the 
"triple  deck  decomposition"  [Hussain  and  Reynolds  (1970)].  This 
means  that  any  instantaneous  quantity  <t>  is  expressed  as 

$  =I(x,y)  +  $p(x,y,ep)  +  <t>(x,y,t)  (2.6) 

where  <i>  is  the  time-mean  value,  i>p  the  periodic  deterministic 
component,  $  the  turbulent  fluctuation  and  the  phase  position 
within  a  cycle.  The  first  step,  in  practice,  is  to  obtain  the 


ensemble-averaged  value  <$>(x,y,9p)  by  averaging  the  instantaneous 

values  at  identical  9  positions  over  a  large  number  of  cycles.  The 

P 

time-mean  value  is  then  obtained  by  averaging  <$>  over  the  complete 

cycle  (0  <.  0p  <  2tt ) .  The  periodic  component  <i>p  is  the  difference 

<$>  _  J.  Finally,  the  turbulent  fluctuation  <j,  is  simply  $  -  <$>. 

This  is  a  statistical  quantity,  for  which  the  phase-averaged  and 

2  ~2 

time-averaged  mean  squared  values,  <$  >  and  $  ,  respectively,  are 
also  calculated  using  the  above  procedure.  For  more  details  of  the 
whole  procedure,  see  Ramaprian  and  Tu  (1982). 


2.5.  Experimental  Conditions 

The  nominal  free-stream  conditions  are  expressed  by  Eq.  (2.1), 
with  U0  =  90  cm/sec  and  e  =  0.40,  and  two  different  frequencies,  f  = 
0.5  Hz  and  f  =  2  Hz.  They  were  attainable  within  +0.5  cm/s  for  the 
velocity,  +  0.0015  for  the  relative  amplitude,  and  ±  0.005  Hz  for  the 
frequency. 

The  turbulent  bursting  frequency  f^  in  the  boundary  layer,  which 
can  be  regarded  as  a  characteristic  frequency  of  turbulence,  was 
estimated,  from  the  criterion  of  Rao  et  al.  (1971)  [see  Eq.  (4.9)]  to 
vary  from  about  13  to  4.5  Hz  from  the  first  to  the  last  measuring 


station. 


It  is  seen  that  the  higher  of  the  two  oscillation 


frequencies  is  comparable  to  the  bursting  frequency,  at  least  for  the 
downstream  stations.  The  fundamental  resonant  acoustic  frequency  of 
the  tunnel  is  estimated  to  be  about  40  Hz.  The  turbulent  energy 
spectrum  was  found  to  contain  only  a  very  small  amount  of  the  energy 
beyond  40  Hz.  Hence,  no  significant  acoustic  interaction  with  the 
turbulence  is  expected. 


The  number  of  sampled  cycles  was  chosen  equal  to  1000  for  the  2- 


Hz  experiment.  In  order  to  keep  the  experimental  time  within 
reasonable  limits,  this  was  reduced  to  250  for  the  0.5-Hz 
experiment.  However,  to  obtain  statistically  stable  averages  for  the 
turbulent  quantities,  the  results  for  each  four  consecutive  phase- 
positions  in  a  cycle  were  averaged,  assigning  the  result  to  the 
intermediate  value  of  the  phase.  This  brought  the  effective  number 
of  samples  again  to  1000. 

Estimations  of  the  (random)  measurement  errors  are  presented  in 
Table  2.1.  They  were  obtained  from  the  scatter  in  the  data  from 
repeated  trials.  Systematic  errors  may  also  occur.  In  the  LDA 
measurements,  the  main  source  of  systematic  error  is  the  deter¬ 
mination  of  the  distance  between  laser  beams  as  they  exit  from  the 
transmitting  optics.  This  quantity  is  used  in  determining  the 
calibration  constant,  and  may  produce  a  deviation  of  the  order  of 
3%.  In  the  HFG  measurements  systematic  errors  are  considered  to  be 
negligible  in  comparison  to  random  ones.  The  estimated  error  in  the 
distance  of  the  measuring  volume  from  the  wall  is  0.4  mm  in  the  outer 
region  (y  >  1  cm)  and  0.12  mm  in  the  inner  one  (1  mm  <  y  _<  1  cm). 
The  (x)  location  of  the  measurement  station  along  the  tunnel,  is 
specified  with  a  tolerance  of  1  cm. 

2.6.  Flow  Quality 

To  assess  the  quality  of  the  flow  in  the  test  section,  the 
following  tests  were  performed: 
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(i)  Check  for  uniformity  and  two-dimensionality:  this  amounts 
to  verifying  the  absence  of  organized  structures  (such  as  vortices) 
in  the  core  of  the  flow  and  determining  the  range  of  influence  of  the 
side  walls. 

(ii)  Check  for  design  conditions:  this  included  verification 
of  Eq.  (2.1),  i.e.  establishing  the  variation  of  the  time-mean  free- 
stream  velocity  and  the  amplitude  and  phase  of  its  oscillatory 
component  along  the  test  section. 

(iii)  Check  for  performance  of  the  tunnel  as  a  steady  flow 

facility:  this  required  comparison  of  the  bounda ry -1 ayer  structure 

to  standard  data  available  in  the  literature. 

The  results  of  each  of  these  tests  are  presented  and  discussed 
in  the  three  following  subsections. 

2.6.1.  Uniformity  and  Two-Dimensionality 

The  uniformity  of  the  flow  in  the  core  region  (beyond  the 
boundary  layers  on  the  walls)  was  verified  in  the  steady  flow  at 
extreme  upstream  and  downstream  stations,  x  =  45  and  190  cm, 
respectively.  The  results  for  the  latter  station  are  presented  in 
Fig.  2.5,  and  are  seen  to  be  very  sati sfactory .  In  fact,  the 
standard  deviation  of  the  velocity  is  only  0.61%  of  the  mean  velocity 
at  this  station,  and  0.64%  at  x  =  45  cm. 

To  check  the  two-dimensionality  of  the  flow  in  the  core  region, 
similar  measurements  were  made,  except  that  they  were  spaced  much 
closer  together.  The  results  for  the  station  x  =-  190  cm  are  shown  in 
Fig.  2.6.  It  is  observed  that  the  flow  is  free  from  any  kind  of 


organized  structure. 


Two-dimensionality  checks  were  also  performed  in  the  boundary 
layer.  This  time,  however,  complete  sets  of  velocity  profiles  were 
measured  at  eight  spanwise  locations  in  the  unsteady  flow  at  the 
higher  frequency  (2  Hz)  at  the  most  downstream  station  (x  =  203 
cm).  The  results  for  the  time-mean  value  of  the  longitudinal 

velocity  are  presented  in  Figs.  2.7  and  2.8  for  the  "left"  (z  <  25 
cm)  and  "right"  side  of  the  test  section  (z  >  25  cm)  respectively, 
when  facing  the  flow.  In  these  figures,  the  parameter  Re  stands 
for  yUg/v-  The  bottommost  distribution  in  both  figures  corresponds 
to  the  plane  z  =  18.9  cm  where  all  the  subsequent  detailed  studies 
were  performed.  It  is  superimposed  as  a  line  on  the  other 
distributions  so  as  to  assess  the  extent  of  two-dimensionality.  It 
is  observed  that  the  agreement  is  good,  at  least  in  the  range  13.6  cm 
<  z  <  29.6  cm.  It  is  necessary  to  mention  that  the  measurements  at 
the  far-right  locations  (largest  z  values)  are  of  lower  quality  due 
to  optical  constraints  imposed  by  the  larger  distance  from  the 
transmitting  optics  and  of  the  LDA.  Furthermore,  a  misalignment  of 
the  right  wall  was  discovered  and  corrected  after  these  measurements 
had  been  completed.  Hence,  an  improvement  in  the  performance  of  the 
flow  on  the  right  side  may  be  expected,  producing  a  more  symmetrical 
behavior.  Similar  results  were  Cutained  for  the  amplitude  and  phase 
of  the  first  harmonic  of  the  longitudinal  velocity,  as  shown  in  Figs. 

2.9  to  2.12.  Results  for  the  time-mean  value  and  amplitude  of 
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oscillation  of  the  turbulent  quantities  <u  >  and  -  <uv>  are  presented 
in  Figs.  2.13  to  2.20.  Again,  the  same  conclusions  hold.  The  flow 

in  the  boundary  layer  can,  therefore,  be  considered  to  be  two- 

dimensional  at  least  over  the  central  30%  of  the  span. 


2.6.2.  Design  Conditions 

The  design  conditions  were  checked  in  both  the  steady  and 
unsteady  regimes.  Figure  2.21  presents  the  longitudinal  variation  of 
the  time-mean  velocity,  relative  to  its  value  at  the  first  measured 
station,  denoted  by  the  subscript  *0*.  Some  of  these  measurements 
were  made  before  the  glass  side  walls  were  removed  (see  section 
2.2).  Two  sets  of  data  points  are  shown  for  the  steady-flow  case, 

corresponding  to  experiments  conducted  on  two  different  days.  They 

agree  well  with  each  other.  It  is  observed  that  the  longitudinal 

variation  of  U  is  within  about  t  3.5%  relative  to  the  mean,  which  is 
acceptable.  The  “peak"  at  x  =  90  cm  seems  to  be  associated  with  some 
waviness  in  the  upper  wall.  In  turn,  it  is  seen  that  the  time-mean 
velocity  for  the  unsteady  flow  at  2  Hz  agrees  very  well  with  the 
steady  flow  data.  Some  deviation  from  the  steady -flow  results  is 
observed  in  the  case  of  the  0.5-Hz  experiment.  This  is  due  to  the 
readjustment  made  in  the  false-wall  position  after  its  reinforcement 
(see  section  2.2).  The  maximum  variation  in  U  is  now  about  ±1.5%. 

The  longitudinal  distribution  of  the  relative  amplitude  of 
oscillation  for  the  unsteady  runs  is  shown  in  Fig.  2.22.  The 

variation  is  of  the  same  order  as  for  the  time-mean  velocity.  Figure 
2.23  presents  the  relative  phase,  which  shows  a  maximum  difference  of 
only  about  4  and  1.20  degrees  for  the  2  and  0.5-Hz  experiments, 
respectively. 

The  higher  harmonic  content  of  the  free-stream  oscillation  was 
found  to  be  of  the  order  of  0.3%  relative  to  the  time-mean  value,  and 
thus  negligible. 


2.6.3.  Steady  Boundary  Layer 


Measurements  of  the  velocity  profile  in  steady  flow  were  made  at 
various  free-stream  velocities  and  at  the  five  measuring  stations, 
namely  x  =  48  cm  (station  1),  x  =  69  cm  (2),  x  =  90  cm  (3),  x  =  142 
cm  (4)  and  x  =  203  cm  (5).  From  the  time-mean  velocity  profiles  the 
integral  parameters,  namely  the  displacement  thickness,  5*,  the 
momentum  thickness,  e,  and  the  shape  parameter,  H  =  6^/0,  were 
calculated.  In  addition,  the  data  were  fitted  by  Coles'  "log  +  wake" 
1  aw, 

—  yu 

rt  - [5-5  ’°9  («) 

l.  % 

where  the  function  W(n),  which  characterizes  the  'wake  component',  is 
taken  as 


for  ^  <  1 

6 


for  ^  >  1 

6 


(2.7) 


W(n)  =  2  (3n2  -  2n3)  (2.8) 

The  fitting  was  performed  by  matching  s*  and  e,  and  using  the 
logarithmic  friction  law  [Eq.  (4,23)],  from  which  the  shear  velocity, 
u*,  the  boundary  layer  thickness,  6,  and  the  wake  strength 
parameter,  n,  are  calculated.  Figure  2.24  shows  the  velocity 


profiles  in  the  inner  or  wall  coordinates  at  station  5  and  the 
corresponding  fit  to  the  Coles-profi le,  Eq.  (2.7).  It  is  observed 
that  the  profiles  do,  in  fact,  follow  very  closely  Coles'  law  at  all 


the  free-stream  velocities.  Very  similar  results  were  obtained  for 
the  remaining  stations.  Table  2.2  presents  the  values  of  the 
parameters  characterizing  the  velocity  profile  at  the  five  measuring 
stations,  for  a  free-stream  velocity  of  Ue  =  90  cm/s.  They  will  be 
useful  as  a  reference  when  analyzing  the  unsteady  flow  results. 

Figure  2.25  shows  the  correlation  between  Cf  and  ReQ  according 
to  the  present  measurements,  and  a  line  representing  the  experimental 
results  of  ten  different  investigations  [Coles  (1962)],  including  one 
performed  at  IIHR  [Landweber  and  Siao  (1958)].  The  agreement  is  very 
good,  except,  perhaps,  at  station  1.  This  may  be  due  to  the  fact 
that  the  flow  was  not  yet  fully  turbulent  at  this  station.  Figure 

2.26  presents  the  variation  of  n  with  Re  ,  together  with  the  "normal 
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curve"  given  by  Coles  (1962).  Except  for  the  points  corresponding  to 
station  3,  which  show  excellent  agreement  with  the  "normal  curve", 
the  present  results  are  consistently  lower.  This  may  be  because  of 
the  relatively  high  free-stream  turbulence  intensity  (1  to  1.2%)  in 
the  present  experiments  [Coles  (1962)],  which,  in  turn,  is  evidently 
due  to  the  absence  of  screens  in  the  calming  chamber  (see  section 
2.2). 

The  four  velocity  profiles  at  station  5  are  replotted  in  Fig. 

2.27  in  the  outer  coordinates  and  compared  with  the  results  of 
Klebanoff  (1954)  corresponding  to  Re0  =  7800.  The  agreement  is 
excellent  for  the  largest  free-stream  velocity,  for  which 

Re  =  5840,  and,  as  expected,  gets  poorer  for  lower  velocities  due 
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to  the  incomplete  development  of  the  wake  component. 


The  distributions  of  the  turbulence  properties,  u  ,  v  and  -  uv 
at  station  5,  normalized  by  the  shear  velocity  u*,  at  station  5,  are 
presented  in  Figs.  2.28  to  2.30,  together  with  the  results  of 
Klebanoff.  The  comparisons  are  considered  to  be  satisfactory  at  the 
higher  free-stream  velocities  and  become  poorer  (as  expected)  at  the 
lower  ones. 

The  results  show  that  the  facility  is,  in  general,  able  to 
produce  a  standard  steady  boundary -layer  flow.  Note  that,  in 
addition,  they  provide  the  data  needed  for  the  "in  situ"  calibration 
of  the  HFG,  during  the  measurement  of  the  wall  shear  stress. 
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Figure  2.6.  Check  of  two-dimensionality  in  the  core  region 


Figure  2.8.  Check  of  two-dimensionality  for  U.  o  , - , 

z  =  18.9  cm;  *,  z  =  24.3  cm;  +.  z  =  29.6  cm;  x, 
z  =  34.9  cm;  z  =  40.2  cm. 
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Figure  2.19.  Check  of  two-dimensionality  for  [-<uv>]. 
Symbols  as  in  Fig.  2.7. 


Figure  2.21.  Longitudinal  distribution  of  the  time-mean  free-stream  velocity. 

A,  v,  steady  flow;  □  ,  unsteady  flow  at  0.5  Hz;Q,  unsteady  flow 


CHAPTER  3 


EXPERIMENT AL  RESULTS 


3.1.  General 

In  this  chapter,  the  main  experimental  results  on  unsteady 
boundary  layers  are  presented  and  discussed.  Emphasis  is  placed  on 
the  comparison  among  the  steady,  quasi -steady  and  unsteady  flows.  In 
the  study  of  the  oscillatory  motion,  the  first  harmonic  of  the 
Fourier  decomposition  of  the  flow  quantities  is  analyzed  in  detail, 
and  only  a  relatively  short  discussion  is  included  on  the  behavior  of 
the  complete  periodic  flow. 

It  is  pointed  out  that,  in  the  interest  of  brevity,  only  typical 
results  are  presented  here.  As  already  mentioned,  all  the  processed 
data  for  both  steady  and  unsteady  flows  are  stored  on  digital  tape 
and  are  available  to  any  interested  reader. 

3.2.  Time -Mean  Flow 

The  time-mean  flow  may  be  expected  to  differ  from  the  steady 
flow  at  the  time-mean  Reynolds  number  due  to  the  nonlinearity  of  the 
problem.  This  nonlinear  behavior  arises,  mathematically  speaking, 
from  the  nonlinearity  of  the  Navier-Stokes  equations.  However,  at 
the  amplitude  of  oscillation  studied  in  the  present  experiments, 
nonlinear  effects  are  found  to  be  practically  insignificant,  in 
agreement  with  previous  investigations. 
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3.2.1.  Velocity  and  Wall  Shear  Stress 

The  time-mean  velocity  for  stations  5  and  1  are  shown, 
respectively,  in  Figs.  3.1  and  3.2  for  both  forcing  frequencies, 
together  with  the  corresponding  quasi-steady  (zero-f requency ) 
velocity  distribution  for  the  same  amplitude  of  oscillation.  The 
latter  distributions  were  obtained  for  each  station  by  averaging  over 
the  steady  distributions  correspondi ng  to  the  complete  range  of  free- 
stream  velocities  included  in  the  oscillation  [i . e . ,  from  (1-e)  U0 

to  (1-e )  u  ].  These  steady  velocity  distributions  were  represented 
by  Coles'  profiles  [Eq.  (2.7)],  whose  parameters  were  related  to  the 
free-stream  velocity  by  means  of  relations  u*(Ue),  n(Ue)  and  5 ( L>e ) 
deduced  from  the  present  steady  flow  measurements.  The  quasi -steady 
distributions  so  obtained  were  found  practically  to  coincide  with  the 
corresponding  steady-flow  velocity  profiles  at  the  time-mean  free- 
stream  velocity.  This  confirms  that  nonlinear  effects  on  the 

velocity  are  not  significant.  On  the  other  hand,  the  time-mean 
velocity  distributions  for  both  frequencies  of  oscillation,  while 
being  very  similar  to  each  other,  are  seen  to  differ  from  the  steady- 
flow  distributions.  This  difference  between  the  steady  and  unsteady 
flow  velocities  is  in  disagreement  with  previous  experimental 
results,  and  has  to  be  explained.  A  first  insight  into  the  nature  of 
this  problem  is  provided  by  the  analysis  of  the  behavior  of  the 

parameters  Reg,  H  and  Cf. 

Figures  3.3  and  3.4  present  the  evolution  of  Refi  and  H  with  Rex 
for  the  steady  flow  at  the  time-mean  velocity  and  the  time-mean 
unsteady  flows.  They  show,  as  expected,  the  same  trends  as  the 


velocity  profiles.  The  values  of  the  parameters  correspond ng  to  the 
two  unsteady  experiments  are  practically  equal  to  each  other,  and 
they  differ  from  the  ones  for  the  steady  flow.  This  difference  tends 
to  disappear  when  moving  downstream.  This  seems  to  indicate  that 
there  is  a  'Reynolds  number  effect',  which  vanishes  asymptotically 
for  large  Reynolds  numbers,  i.e.,  when  the  flow  is  completely 
developed.  The  results  for  the  time-mean  wall  shear  stress  are  shown 
in  Fig.  3.5.  It  seems  that  the  time-mean  unsteady  value  for  the 
smaller  frequency  is  slightly  larger  than  for  the  higher  frequency, 
and  very  close  to  the  wall  shear  stress  for  the  steady  flow  at  the 
time-mean  free-stream  velocity.  This  latter  value,  in  turn,  is 
always  smaller  than  the  quasi -steady  value.  Thus,  the  trend  seems  to 
be  a  slight  decrease  in  time-mean  wall  shear  stress  with  the 
frequency.  It  is  interesting  to  note  that  this  is  opposite  to  the 
trend  observed  by  Ramaprian  and  Tu  (1982)  in  pipe  flow.  In  any  case, 
the  magnitude  of  the  decrease  is  too  small  to  be  quantified  from  the 
present  measurements,  which  have  an  uncertainty  of  ±  1%. 

From  the  experimental  results  shown,  it  is  obvious  that  the 
outer  and  inner  layers  are  responding  quite  differently  to  the 
imposed  oscillation.  On  the  one  hand,  the  inner  flow  (whose  behavior 
is  illustrated  by  that  of  Cf)  does  not  seem  to  be  very  much 
affected.  On  the  other  hand,  the  outer  flow  (whose  response  is 

characteri /ed  by  those  of  Re  and  H)  is  s i gni f i cantly  changed  by 

0 

unsteadiness,  though  the  differences  tend  to  disappear  for  large 
ReK.  This  mismatch  is  well  illustrated  in  Fig.  3.6,  where  the 

relation  between  Cf  an,j  for  the  unsteady  fl  ow  experiments  is 
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compared  with  the  experimental  data  reviewed  by  Coles  (1962).  It  is 
seen  that  they  do  not  follow  the  "normal"  law  for  a  standard  steady 
turbulent  boundary  layer,  except  for  the  largest  ReQ  values. 

In  order  to  understand  the  details  of  behavior  of  the  outer  and 
inner  flow  in  the  unsteady  flow  regime,  it  is  necessary  to  analyze 
the  complete  velocity  profile.  Figures  3.7  and  3.8  present  the 
velocity  distributions  in  the  inner  coordinates  for  the  five  stations 
and  the  two  frequencies.  The  measured  value  of  the  time-mean  shear 
velocity  is  used  as  the  scaling  velocity.  Lines  representing  the 
universal  logarithmic  law  (u+  =  5.5  log  yf  +  5.5)  and  the  quasi¬ 
steady  velocity  distribution  are  also  shown.  It  is  observed  that  the 
time-mean  unsteady  velocity  profiles  do  not  follow  the  universal  law, 
except  asymptotically  at  large  Reynolds  numbers  (i.e.,  at  the  more 
downstream  stations).  However,  they  show  a  logarithmic  region,  but 
with  a  slightly  different  slope  and  a  significantly  different 
additive  constant. 

In  fact,  these  time-mean  distributions,  like  their  steady 
counterparts,  also  belong  to  the  Coles  family  of  velocity  profiles 
[see  Eq.  (2.7)].  The  corresponding  parameters  (Cf,  6  and  it),  for 
each  flow  and  each  station,  are  presented  in  Table  3.1,  together  with 
the  ones  for  the  steady  flow  at  the  time-mean  velocity,  as  a 
reference.  Note  that,  in  the  case  of  the  unsteady  flows,  the 
physical  interpretation  of  Coles'  parameters  is  not  obvious  except 
for  (5,  which  can  be  considered  as  the  boundary  layer  thickness  for 
the  time-mean  flow.  It  is  this  quantity  that  is  used  in  Figs.  3.9 
and  3.10  to  present  the  vel oci ty -defect  distribution  in  the  outer 


coordinates.  Also  shown  are  the  quasi -steady  profiles  and  a  line 

representing  Klebanoff's  data  (1954),  considered  to  be  the  standard 
distribution  for  steady  flow  at  large  Reynolds  numbers.  It  is  seen 
that  both  the  unsteady  and  quasi-steady  distributions  approach 

Klebanoff’s  profile  for  large  Rex,  the  unsteady  profiles  doing  so 
faster,  though  nonmonotoni cal ly .  In  any  case,  the  approach  is  from 
'below',  showing  that  there  is  development  of  the  wake  component  in 
the  downstream  direction. 

It  is  interesting  to  go  back  and  study  the  significance  of 
Coles'  parameters  for  the  time-mean  velocity  profile  in  unsteady 
flow.  The  values  of  the  parameter  Cf  are  also  shown  in  Fig.  3.5.  It 
is  observed  that,  except  at  the  last  stations,  they  do  not  correspond 
to  the  correct  value  of  the  time-mean  wall  shear  stress  (unlike  in 
steady  flow).  However,  when  plotted  in  the  Coles-plot  of 

vs.  CfRe0  (see  Fig.  3.6),  they  fall  on  the  'normal'  line.  Hence, 

they  can  be  interpreted  as  the  skin  friction  coefficient  for  a 
'normal'  steady  turbulent  boundary  layer  at  the  given  Re0.  As 
for  n,  it  can  still  be  regarded  as  a  measure  of  the  strength  of  the 
wake  component  relative  to  the  logarithmic  region.  Its  variation 
with  Re0  is  presented  in  Fig.  3.11.  Comparison  with  Fig.  2.26  shows 
that  the  distribution  is  'normal'  by  the  present  standards. 

These  last  results  demonstrate  that,  in  the  unsteady  boundary 
layer,  the  mean  outer  flow  can  be  considered  to  be  the  same  as  that 
in  a  'normal'  steady  boundary  layer  at  the  same  Re  .  From  this  point 
of  view,  the  observed  difference  between  the  velocity  distributions 
at  a  given  location,  in  steady  and  unsteady  flows,  can  be  interpreted 


as  being  due  to  a  relative  shift  in  the  virtual  origins  for  the  two 
cases.  This  is  also  strongly  suggested  by  Fig.  3.3,  from  which  one 
can  obtain,  by  extrapolation,  the  approximate  position,  xv,  of  the 
virtual  origin.  It  is  found  that  xv  -  7  cm  for  the  steady  flow  and 
xv  =  -  75  cm  for  both  the  time-mean  unsteady  flows.  Hence,  the 
effect  of  the  imposed  oscillation  on  the  outer  flow  is  to  promote  a 
rapid  development  towards  its  asymptotic  structure.  Note  that  the 
response  seems  to  be  quantitatively  similar  for  both  frequencies. 
However,  from  Table  3.1  or  Fig.  3.11  it  is  observed  that  n  actually 
varies  significantly  from  one  frequency  to  the  other  at  a  given 
station.  Since  both  the  unsteady  profiles  have  similar  values 

of  Re.  at  all  stations,  this  difference  must  he  due  to  the  effect  of 

0 

unsteadiness.  We  shall  call  it  the  "frequency  effect".  It  is 
plausible  (and  consistent  with  other  results,  as  shown  in  the  next 
section)  to  assume  that  frequency  effects  manifest  themselves  as  an 
increase  in  n.  Hence,  they  seem  to  be  important  at  the  last  station 
for  f  =  0.5  Hz  and  at  the  first  three  stations  for  f  =  2  Hz. 

In  conclusion,  it  can  be  said  that  the  imposed  oscillation  acts 
very  differently  on  the  inner  and  outer  layers.  The  inner  flow 


responds 

rather  weakly.  The 

outer 

fl  ow, 

on  the 

other  hand,  is 

strongly 

affected.  The  changes 

in  the  outer 

flow  can 

be  interpreted. 

in  general,  as  a  combination 

of 

rather 

strong 

Reynolds  number 

effects. 

which  are  due  to  shift 

of  the 

vi rtual 

origin  of  the 

turbulent  boundary  layer,  and  relatively  mild  frequency  effects  which 
increase  the  wake-component  strength.  The  former  are  almost 
independent  of  frequency  and  decrease  in  the  downstream  direction, 
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i.e.,  with  increasing  length  Reynolds  numbers.  It  is  seen  that  the 
main  differences,  between  the  time-mean  velocity  distribution  and 
that  for  steady  flow  at  the  time-mean  free-stream  velocity  observed 
in  the  present  experiments,  is  due  to  their  different  'history1.  In 
the  limit  of  large  Reynolds  number,  they  nearly  coincide  with  each 
other  (in  agreement  with  previous  experimental  results),  except  for 
small  residual  effects  which  depend  on  the  oscillatory  frequency. 
Some  correlation  (though  not  an  explanation)  will  be  found  for  these 
frequency  effects  later  in  this  work  (see  chapter  5). 

3.2.2.  Turbulence  Properties 

The  distributions  of  the  mean-square  longitudinal  turbulent 
~2 

intensity,  u  ,  for  the  five  stations  and  for  both  the  frequencies, 
are  presented  in  Figs.  3.12  to  3.16  in  the  outer  coordinates.  Also 
shown  are  the  corresponding  distributions  in  steady  flow  at  the  time- 
mean  velocity,  and  the  data  of  Klebanoff  (1954),  considered  to  be  the 
asymptotic  limit  for  large  Reynolds  numbers.  It  is,  in  fact,  seen 
that  this  limit  is  slowly  but  continuously  approached  by  the  present 
steady-flow  data  for  increasing  values  of  Rex. 

The  behavior  of  the  distributions  in  the  unsteady  flows  is 
rather  interesting  when  compared  with  the  steady  flow.  In  fact, 
trends  not  previously  observed  can  be  clearly  distinguished.  First, 
it  is  to  be  noted  that,  at  large  Reynolds  numbers,  the  steady-flow 
profiles  are  not  very  sensitive  to  Reynolds -number  variations.  This 
is  related  to  the  above-mentioned  insensitivity  f  >  nonlinear  effects, 
and  can  be  checked  by  comparing  the  distributions  for  stations  3,  4, 


« 
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and  5,  which  are  very  similar  to  each  other.  Significant  changes 
occur,  however,  for  stations  1  and  2,  where  Re ^  is  relatively 
small.  Hence,  for  these  stations  an  additional  line  is  drawn, 
corresponding  to  the  distribution  in  steady  flow  at  station  3.  This 
may  be  used  as  the  reference  steady-flow  distribution  for  a  value 
of  Reg  similar  to  those  for  the  unsteady -fl ow  profiles  at  those 
stations.  The  difference  between  this  reference  distribution  and  the 
actual  steady-flow  profile  is  a  Reynolds  number  effect.  In  turn,  any 
departure  of  the  unsteady -fl ow  distribution  from  the  reference 
profile  has  to  be  considered  as  a  frequency  effect.  It  is  observed, 
from  the  figures,  that  frequency  effects  show  up  as  a  positive 
departure,  at  the  two  downstream  stations  for  f  =  0.5  Hz  and  at  the 
two  upstream  stations  for  f  =  2  Hz.  These  differences,  though  not 
very  large,  are  significant  and  consistent  with  the  discussion  on  the 
effects  of  unsteadiness  on  the  velocity  profile. 

Similar  trends  are  observed  for  the  mean-square  transverse 
turbulent  intensity  v  .  Figure  3.17  presents  the  results  for  station 
5,  where  frequency  effects  are  seen  to  be  present  for  f  =  0.5  Hz, 
while  Fig.  3.18  shows  the  distributions  for  station  3,  where  very 
significant  departures  due  to  frequency  effects  are  observed  for  f  = 

2  Hz.  Notice  the  tendency  of  the  v  -distributions  to  increase  very 
close  to  the  wall.  This  is  believed  to  be  associated  with  the 
(small)  oscillation  of  the  test  surface.  Figures  3.19  and  3.20 
illustrate  the  behavior  of  the  turbulent  Reynolds  shear  stress.  Note 
that,  now,  no  significant  frequency  effects  are  observed  for  station 

3  at  any  frequency. 


3. 3.  Osci 1 1 atory  Motion 

Only  the  behavior  of  the  amplitude  and  phase  of  the  first 
harmonic  of  the  flow  quantities  is  studied  in  the  following  sections, 
as  it  is  the  most  energetic  component  of  the  oscillatory  motion.  A 
brief  discussion  of  the  complete  ensembl e-averaged  flow  and  its 
departure  from  a  sinusoidal  response  is  presented  in  a  later  section. 

3.3.1.  Velocity 

Stations  5  and  3  are  selected  as  representative  of  the  behavior 
of  the  oscillatory  motion.  The  correspondi ny  results  for  the 
amplitude  of  oscillation  of  the  velocity  are  shown,  respectively,  in 
Figs.  3.21  and  3.22.  An  overshoot,  typical  of  periodic  flow,  can  be 
seen  in  each  case.  Its  magnitude  is  of  the  order  of  7-8%  of  the 
free-stream  amplitude,  and  decreases  slowly  with  the  position  along 
the  plate  and  the  oscillation  frequency.  The  extent  of  amplitude 
variation,  however,  varies  significantly  both  with  x  and  f.  Beyond 
the  amplitude  overshoot,  and  closer  to  the  outer  edge  of  the  boundary 
layer,  the  flow  undergoes  constant  amplitude  oscillations  at  the 
free-stream  amplitude,  i.e.,  it  is  an  inviscid  or  slug-flow 
oscillation.  The  unsteady  viscous  effects  can  be  considered  to  be 
confined  within  the  zone  of  significant  amplitude  variation.  A 
measure  of  the  extent  of  unsteady  viscous  effects  is  provided  by  the 
distance,  y^,  from  the  wall  to  the  point  of  maximum  amplitude.  This 
distance,  referred  to  the  local  boundary  layer  thickness,  is 
presented  in  Table  3.2  for  each  station  and  the  two  frequencies. 
Also  the  corresponding  relative  btokes-layer  thickness,  namely 
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/  — /S,  is  tabulated  as  a  reference,  which  measures  the  extent  of 
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unsteady  viscous  effects  in  a  hypothetical  periodic  laminar  boundary 
layer  of  the  same  thickness  and  oscillated  at  the  same  frequency.  It 
is  observed  that  yM/6  decreases  in  the  downstream  direction,  for  a 
given  frequency,  and  also  decreases  with  frequency,  for  a  given 
station.  However,  it  is  interesting  to  note  that  when  y^  is 
normalized  by  u*/&j,  it  remains  essentially  constant  with  both 
frequency  and  position  along  the  plate,  as  shown  in  the  last  column 
of  Table  3.2.  In  any  case,  y^/6  is  always  very  much  larger  than  the 
relative  Stokes-layer  thickness,  showing  that  unsteady  viscous 
effects  propagate  over  a  much  wider  region  in  turbulent  flow.  These 
results  are  in  agreement  with  previous  investigations.  The  table 
also  suggests  that,  from  the  point  of  view  of  unsteady  viscous 
effects,  there  is  a  continuous  transition  between  what  happens  at 
station  5  at  the  lower  frequency  and  at  station  1  at  the  higher 
one.  All  these  trends  will  be  shown  to  be  consistent  with  the  theory 
developed  in  the  next  chapter. 

Figures  3.23  and  3.24  present  the  distribution  of  the  phase 
lead,  relative  to  the  free-stream  velocity,  of  the  velocity 
oscillation.  It  increases  monotoni cal  ly  towards  the  wall,  in  some 
cases  after  attaining  a  small  phase  lag.  This  agrees  with  previous 
experimental  results.  At  the  closest  point  (at  about  1  mm  from  the 
wall)  the  phase  lead  is,  for  all  stations,  in  the  range  10-13  degrees 
for  f  =  0.5  Hz  and  7-9  degrees  for  f  -  2  Hz.  The  significant  phase 
variations  extend  over,  roughly,  the  same  distance  as  the  amplitude 
variations.  Also  shown  in  these  figures  are  the  phases  of  the 


corresponding  wall  shear  stress  (see  below),  except  for  station  5  at 
f  =  2  Hz,  in  which  case  it  goes  out  of  the  range  of  this  figure. 
They  suggest  that  some  nonmonotonic  behavior  may  occur  close  to  the 
wall  for  the  1  ow-f requency  cases.  However,  measurements  could  not  be 
obtained  closer  than  1  mm  from  the  wall  in  the  present  experiments. 

3.3.2.  Turbulence  Properties 

The  amplitude  of  oscillation  of  the  turbulence  properties  <u^>, 

p  — 

<v  >  and  -<uv>  are  normalized  by  2[<Ue>]  U  .  The  logic  behind  this 

nondimensional  i  zation  is  that  the  resulting  quantity,  say  the 

2 

amplitude  of  <u  >,  should  be  essentially  equal,  in  quasi-steady  flow, 

~~2  —? 

to  the  steady-flow  value  of  u  /Ue  (see  chapter  4).  Hence,  any 

departure  from  the  steady-flow  distribution  indicates  that  the 

p 

quantity  <u  >  is  not  varying  in  a  quasi-steady  manner.  Figures  3.25 
and  3.26  present  the  amplitude  of  oscillation  for  <u2>.  Departures 
from  the  quasi -steady  behavior  are,  as  expected,  very  important  for 

both  the  frequencies  (compare  with  Figs.  3.12  and  3.14).  On  the 

other  hand,  the  region  of  significant  variation  of  the  amplitude  are 
now  more  extended  (roughly  doubled)  than  for  the  velocity.  Beyond 
this  region  the  amplitude  is  essentially  zero,  indicating  a  state  of 
frozen  turbulence,  in  agreement  with  previous  investigations.  For 
station  3  (Fig.  3.26)  and  f  -  0.5  Hz,  a  secondary  peak  in  the 
amplitude  is  observed  close  to  the  outer  edge  of  the  boundary 
layer.  The  origin  of  this  disturbance,  which  also  shows  up  at  other 
locations  and  at  the  higher  frequency,  could  not  be  sat i sfactori ly 
expl ai ned. 
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The  phase  of  <u  >,  relative  to  the  free-stream  velocity,  is 
presented  in  Figs.  3.27  and  3.28.  It  is  actually  interpreted  as  a 
phase  lag,  which  tends  to  zero  at  or  close  to  the  wall  (where  the 
turbulence  is  produced),  and  increases  to  360  degrees  at  the  outer 
edge  of  the  region  with  strong  unsteady  viscous  effects,  which  will 
be  called  the  "unsteady  layer"  (indicated  by  a  vertical  line  for  f  = 
2  Hz).  The  scatter  observed  beyond  this,  for  f  =  2  Hz,  is  a  result 
of  the  very  small  (spurious)  amplitude  of  oscillation  in  that  region. 

Similar  results  are  obtained  for  <v^>,  as  depicted  in  Figs.  3.29 
to  3.32.  This  time,  however,  the  unsteady  viscous  effects  extend 
over  a  slightly  smaller  region.  Note  the  anomaly  in  the  phase 
distributions  at  f  =  2  Hz  over  the  outer  region  (Figs.  3.31  and 
3.32).  The  amplitude  of  this  disturbance,  however,  is  very  small,  as 
seen  from  Figs.  3.29  and  3.30. 

Lastly,  Figs.  3.33  to  3.36  present  the  amplitude  and  phase  of 
the  Reynolds  shear  stress,  which  show  a  behavior  very  much  analogous 

O 

to  that  of  <u  >. 

The  amplitude  and  phase  distributions  of  the  turbulence 
properties  are  in  general  agreement  with  those  measured  by  Ramaprian 
and  Tu  (1982)  and  Mizushina  et  al.,  (1975)  in  periodic  pipe  flow.  It 
is  seen  that,  as  observed  by  these  investigators,  the  phase  lag  of 
<u^>  increases  approximately  linearly  through  the  unsteady  layer. 
This  can  be  interpreted  as  if  the  disturbance  in  the  turbulence 
(generated  near  the  wall)  imposed  by  the  oscil’ation  diffuses  outward 
with  a  constant  velocity.  Now,  the  distance  travelled  by  the 
disturbance  in  one  oscillation  period,  due  to  diffusion,  corresponds 


to  the  thickness  of  the  unsteady  layer.  If  this  thickness  can  be 
regarded  to  be  characteri zed  by  yM,  the  relation  <4yM/u*  ~  1  indicated 
by  Table  3.2  suggests  that  the  diffusion  velocity  is  constant  and  is 
roughly  equal  to  u*.  Another  significant  point  to  be  noted  is  that 
both  <u 2>  and  -  <uv>  are  out  of  phase  with  <U>  over  most  part  of  the 
unsteady  layer,  a  result  which  is  in  agreement  with  the  observations 
of  Ramaprian  and  Tu  (1982).  This  points  to  the  limitation  of  quasi¬ 
steady  eddy  viscosity  models  that  relate  <uv>  to  <U>.  A  careful 
examination  of  Figures  3.27,  3.28,  3.35  and  3.36  also  shows  that  <u^> 
and  -<uv>  have  a  relative  phase  difference  which  varies  across  the 
unsteady  layer.  This  is,  again,  in  agreement  with  the  earlier 
results  in  pipe  flow  and  shows  that  the  turbulence  structure  is  not 
in  'equilibrium'  during  the  oscillation  cycle. 

3.3.3.  Wal 1  Shear  Stress 

Figure  3.37  presents  the  amplitude  of  oscillation  of  the  wall 
shear  stress,  normalized  in  the  same  manner  as  the  turbulence 
properties.  Comparison  with  the  time-mean  results,  represented  by 
the  shaded  region,  shows  that  there  is  some  departure  from  the  quasi - 
steady  distribution  for  the  amplitude  at  all  the  stations.  This 
departure  is  larger,  as  expected,  for  the  higher  frequency.  However, 
the  deviation  from  quasi-steady  flow  is  not  as  strong  as  that 
observed  in  the  case  of  the  turbulence  properties  in  the  outer  layer. 

The  phase  of  the  wall  shear  stress  is  the  more  reliable  quantity 
measured  with  the  MFC.  The  reason  is  that  it  is  not  very  sensitive 
to  the  calibration  constants.  Figure  3.38  presents  the  results  for 
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the  four  measuring  stations.  It  is  observed  that  the  phase  lead 

increases  in  the  downstream  direction  for  a  given  frequency,  and  also 
increases  with  the  frequency  for  a  given  station. 

Additional  experiments  were  done  with  the  HFG  at  station  5 
varying  the  frequency  of  oscillation,  but  using  the  sleeve  designed 
for  f  =  0.5  Hz  (see  section  2.2).  The  results  for  the  time-mean  wall 
shear  stress  are  shown  in  Fig.  3.39.  A  weak  tendency  to  decrease 
with  the  frequency  is  observed.  The  quasi-steady  value  corresponding 
to  this  station  is  also  shown  in  this  figure,  and  seems  to  agree  well 
with  the  general  trend.  The  amplitude  of  oscillation  is  shown  in 

Fig.  3.40,  where  also  the  time-mean  results  are  represented  as  a 
shaded  region.  Again,  there  is  a  weak  departure  from  a  quasi-steady 
amplitude  distribution.  The  phase  variation  is  presented  in  Fig. 

3.41,  where  results  obtained  using  the  sleeve  designed  for  f  =  2  Hz 
are  also  shown.  A  monotonic  and  roughly  linear  increase  with  the 

frequency  is  observed. 

3.4.  Ensemble-Averaged  Flow 

Though  the  harmonic  motion  with  the  forcing  frequency  is  the 
most  energetic  component,  higher  harmonic  motions  are  not  negligible 
for  the  present  amplitude  of  oscillation.  Only  the  ensemble  averaged 
velocity  can  be  considered  sinusoidal  with  little  error,  as 
illustrated  in  Fig.  3.42,  where  the  phasewise  variation  of  the  free- 
stream  velocity  is  shown.  For  any  location  across  the  boundary 
layer,  higher  harmonic  contributions  to  <U>  were  found  to  be  less 
than  1%  of  the  fundamental. 
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This  is  not  the  case  for  the  turbulence  properties.  Figures 
3.43  to  3.45  present  their  variation  with  phase  position  for 
different  locations  relative  to  the  wall,  at  station  5  and  f  =  2 
Hz.  One  location  (y/6  =  .47)  corresponds  to  the  region  where  the 
turbulence  structure  is  'frozen';  another  (y/5  .068)  to  the  zone  of 
maximum  oscillation  amplitude;  and  a  third  (y/5  =  .023)  to  tne 
measurement  point  closest  to  the  wall.  The  contribution  from  the 
higher  harmonics,  relative  to  the  fundamental,  is  of  the  order  of  15- 
20%  at  y/6  =  .07  and  increases  up  to  about  40%  at  y/6  =  .02.  Similar 
results  were  observed  in  the  lower  frequency  experiment.  Higher 
harmonic  contributions  in  this  case  ranged  from  10-20%  at  the  zone  of 
maximum  amplitude  (y/6  =  .22)  to  about  30%  at  the  closest  measuring 
point. 

Figure  3.46  shows  the  phasewise  variation  of  the  ensemble 
averaged  wall  shear  stress  in  some  typical  cases.  Higher  harmonics 

in  <t  >  were  found  to  be  of  the  order  of  7-10%  for  f  =  0.5  Hz  and  15- 

w 

20%  for  f  =  2Hz,  relative  to  the  fundamental. 


3.5.  Derived  Quantities 

The  experimental  information  can  be  used  to  calculate  some 
secondary  quantities  that  are  important  from  a  theoretical  and/or 
practical  point  of  view.  Some  of  these  are  presented  and  briefly 
discussed  next. 

Figures  3.47  to  3.49  show  the  phasewise  variation  of  the 


integral  parameters  for  a  few  phase  positions  along  the  cycle,  at 
station  5  and  for  both  frequencies.  It  is  seen  that  the  time-mean 
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value  as  well  as  the  amplitude  of  oscillation  are,  for  each 

parameter,  very  similar  for  both  the  frequencies.  The  relative 
amplitude  is  roughly  43%  for  the  displacement  thickness,  33%  for  the 
momentum  thickness  and  only  12%  for  the  shape  parameter.  In  quasi - 
steady  flow,  these  integral  parameters  are  180  degrees  out  of  phase 
with  the  free-stream  velocity.  It  is  observed  that  this  is  still 
roughly  the  case  for  the  0.5-Hz  experiment,  but  that  for  the  2-Hz 
experiment  there  seems  to  be  a  phase  lead,  relative  to  the  quasi - 
steady  case,  of  about  30  degrees.  It  is  seen  that  there  is 
significant  distortion  from  the  sinusoidal  behavior. 

It  was  found,  from  a  detailed  study  of  the  data,  that  the 

ensemble-averaged  velocity  profiles  for  given  phase  positions  do  not 
belong  to  the  Coles'  family.  However,  it  seems  to  be  still  possible 
to  distinguish  a  logarithmic  region  close  to  the  wall  which  matches 
with  some  of  the  lines  in  a  Cl auser-pl ot ,  though  smetimes  this  is  a 
difficult  and  ambiguous  task.  The  "skin  friction  coefficient"  so 
obtained,  for  station  5,  is  presented  in  Fig.  3.50.  Compared  to  the 
results  of  the  wall  shear  stress  measurements ,  their  mean  values  and 

amplitudes  are  low,  specially  for  the  0.5  Hz  case.  This  means  that 

the  conventional  Clauser-plot  cannot  be  used  to  obtain  the  phasewise 
distribution  of  the  wall  shear  stress. 

Another  important  and  useful  quantity  is  the  eddy  viscosity. 
Its  ensembl e-averaged  value  is  defined  as 

,  v  -  -<uv> 

<vt>  ■=  ~3  <11  > 

ay 


(3.1) 


To  obtain  a  meaningful  estimate  of  this  value,  it  is  crucial  10 
smooth  the  variation  of  <l>  with  y  before  proceeding  with  the 
differentiation.  This  is  difficult  and  very  laborious  to  do  with  the 
whole  set  of  ensemble-averaged  velocity  profiles.  A  simpler  method 
can  be  used  by  invoking  the  smal 1 -ampl i tude  approximation.  Utilizing 
the  complex  notation,  one  can  write 


<U>  =  U  +■  cl^e’^  (3.2) 

-<uv>  =  r  =  t  +  sT^e'^  (3.3) 

<'vt'>  =  vt+£vtle^lil  (3.4) 

Introducing  Eqs.  (3.2)  to  (3.4)  into  Eq.  (3.1),  the  following 
expressions  for  the  time-mean  value  and  the  complex  amplitude  of  the 
eddy  viscosity  follow: 


V 


t 


3lJ 

ay 


(3.5) 


utl 


-  "“1 
1 l"ut  ay 

ay 


(3.6) 


Thus,  the  smoothing  procedure  is  now  only  needed  for  U  and  (which 
includes  both  the  "in-phase"  ana  "out -of -phase"  components).  Coles' 
c )  t  is  used  for  tt  while  Ui  is  represented  by  a  piecewise  pol  ynomial 


fit.  Figure  3.51  presents  the  time-mean  distribution  for  station  5 
and  both  frequencies.  A  line  is  also  shown,  representing  the 
distribution  of  in  steady  flow.  As  in  the  case  of  the  turbulent 

shear  stress,  the  time-mean  distribution  of  for  the  higher 
frequency  agrees  very  well  with  that  in  steady  flow,  but  the 
distribution  for  f  =  0.5  Hz  shows  a  departure  due  to  frequency 
effects.  Figures  3.52  and  3.53  present  the  amplitude  and  phase  of 
the  oscillatory  component.  Note,  that  for  f  -  2  Hz  the  eddy 
viscosity  is  frozen  in  the  outer  region  (as  is  to  be  expected  since 
both  <-■>  and  3<U>/3y  are  frozen). 


Table  3.1 


Coles*  Parameters  for  the  Time-Mean  Unsteady  Flows 
and  the  Steady  Flow  at  the  Time-Mean  Free-Stream  Velocity 
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Figure  3.6.  Correlation  between  Cf  and  Re„  For  unsteady  flow 
□,  f  =  0.5  Hz;  0,  f  =  2  Hz;  +,  Coles  parameter 
for  f  =  0.5  Hz;  *,  Coles  parameter  for  f  =  2  Hz; 
- ,  'normal  curve'  [Coles  (1962)]. 
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Figure  3.9.  Time-mean  velocity  profiles  in  the  outer 

coordinates  for  f  =  0.5  Hz;  - ,  Klebanoff 

data  (1954).  Other  symbols  as  in  Fig,  3.7. 
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Time-mean  distribution  for  <u  at  station  5 
*,  f  =  0.5  Hz;  0  ,  f  =  2  Hz;  ,  steady  flow 
- ,  Klebanoff's  data  (1954). 
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19.  Time-mean  distribution  for  -  <uv>  at  station  5 
Symbols  as  in  Fig.  3.12. 


Figure  3.21.  Amplitude  distribution  for  the  velocity  at 
station  5.  *,  f=0.5Hz;o,f=2Hz. 


Phase  distribution  for  <u  >  at  station  3. 

edge  of  unsteady  layer.  Other  symbols 
as  in  Fig.  3.21 . 


Phase  distribution  for  -  <uv>  at  station  5. 
— ,  edge  of  unsteady  layer.  Other  symbols 
as  in  Fig.  3.21. 


Frequency  variation  of  Cf  at  station  5.  A,  quasi-steady  value 


Phasewise  variation  of  [<Cf>]  at  f  =  2  Hz. 


0.  80.  160.  270.  360. 

Gp  (deep 

Figure  3.47.  Phasewise  distribution  of  <o*>.  0,  f  =  0.5  Hz; 

0,  f  =  2  Hz. 
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Figure  3.50.  Phasewise  distribution  of <C^>.  0,  f  =  2  Hz,  from 

Clauser-plot;  •  ,  f  =  0.5  Hz,Tfrom  Clauser-plot;  *, 
f  =  2  Hz,  measured;  +,  f  =  0.5  Hz,  measured. 


Figure  3.53.  Phase  distribution  for  <v 
layer.  Other  symbols  as 
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PERIODIC  T  i L F NT  FLOWS 
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4.1.  General 

The  stumbling  block  in  the  a": '/sis  of  turbulent  flows,  both  in 
the  steady  and  unsteady  regimes,  is  the  well-known  closure  problem. 
This  arises  as  a  consequence  of  the  statistical  approach,  that 
provides  a  hierarchy  of  equations  in  which  correlations  of  a  given 
order  are  related  to  the  next-order  correlations.  As  an  attempt  to 
overcome  this  difficulty,  the  so-called  turbulence  models  have  been 
introduced.  Essentially,  they  amount  to  truncating  the  sequence  of 
equations  at  a  certain  level,  and  relating  the  remaining  highest- 
order  correlations  to  the  lower-order  ones  through  "models"  based  on 
ad  hoc  physical  arguments  and  experimental  data  [see  Bradshaw 
(1976)]. 

Though  useful  for  engineering  purposes,  the  turbulence  models  do 
not  seem  to  be  the  appropriate  tool  to  provide  new  physical  insight 
into  the  turbulence  process.  Such  insight  has  been  obtained 

primarily  from  experiments  and,  to  some  extent,  from  simple,  yet 
powerful,  dimensional  analysis.  In  fact,  a  combination  of  a  minimum 
of  experimental  information,  dimensional  analysis  and  singular 
perturbation  theory  has  proved  useful,  in  recent  times,  in  providing 
quite  general  asymptotic  results  for  large  Reynol  is  numbers,  in  the 
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form  of  similarity  laws  for  many  categories  of  steady  turbulent  shear 
flows  [see  Tennekes  and  Lumley  (1980)].  The  power  of  this  approach 
(and,  perhaps,  its  major  limitation)  is  that  it  works  on  the  open 
system  of  equations,  thus  avoiding  the  anomalies  introduced  by  a 
specific  turbulence  model.  While  in  the  case  of  steady  flows  such 
analysis  has  merely  reconfirmed  already  well-established  experimental 
knowledge,  it  is  used  in  the  present  work  as  a  predicting  tool  to 
obtain  similarity  laws  for  unsteady  turbulent  flows.  Some  of  these 
results  are  verified  in  the  next  chapter,  using  the  available 
experimental  information. 

In  analyzing  the  behavior  of  the  velocity  in  unsteady  (periodic) 
flows,  the  tendency  has  been  to  search  for  similarity  laws  (e.g.,  law 
of  the  wall,  logarithmic  law,  etc.)  in  the  instantaneous  (ensemble- 
averaged)  velocity  profiles  [Cousteix  et  al.  (1931)],  with  little 
success.  In  the  present  work,  instead,  the  velocity  profile  is 
decomposed  into  its  time -mean  and  oscillatory  components,  and  each 
one  of  them  is  studied  separately. 


4.2.  Equations 


The  analysis  is  started  from  the  ensemble-averaged  unsteady 
turbulent  boundary-layer  equations  for  an  incompressible,  two- 
dimensional  flow  [Patel  and  Nash  (1971)] 


3<U>  3<V> 

3x  3y 


(4.1) 


!<M>  +  <u>  ®<y>  +  <v>  i<J!> 

3t  3X  3y 


l  Kpj>  32<U>  3 <t> 

p  ax  +  v  +  ~y~ 


(4.2) 
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1  3<p> 
p  3X 


+  U 


e  3x 


(4.3) 


The  inclusion  of  the  unsteady  terms  in  Eqs.  (4.2)  and  (4.3)  does  not 

involve  any  additional  hypothesis  beyond  the  classical  boundary-layer 

assumption  for  steady  turbulent  flows.  Briefly,  this  assumption 

states  that  the  length  scale  of  the  motion  normal  to  the  wall  is  much 

smaller  than  that  along  the  wall.  As  a  consequence,  diffusion  in  the 

longitudinal  direction  can  be  neglected,  and  the  pressure  can  be 

considered  to  remain  constant  in  the  transverse  direction.  Actually, 

in  turbulent  flow,  this  last  conclusion  requires  one  to  neglect  also 

2 

the  contribution  of  the  normal  Reynolds  stress  <v  >,  which  can  be 

done  safely  except  near  the  separation  point.  Then,  the  longitudinal 

pressure  gradient  is  a  known  quantity  in  boundary-layer  theory,  and 

is  related  to  the  free-stream  velocity  by  Eq.  (4.3).  Likewise,  in 

2 

Eq.  (4.2)  the  contribution  due  to  the  normal  turbulent  stress  <u  > 
has  been  neglected,  for  which  the  same  considerations  apply  as  for 
<v^>.  The  above  equations  are  subject  to  the  following  boundary 
conditions  • 


<’)>  =  0,  <V>  =  0  at  y  =  0 

<U>  •*  L)a  as  y  ■*  °°  (A. 4) 

< r >  =  0  at  y  =  0 

The  upstream  boundary  conditions,  i.e.,  the  specification  if  <'J'>  at 
tie  first  station,  is  somewhat  complicated. 


Howovv',  the  analysis 


will  concentrate  on  the  search  for  asymptotic  sel f-simi 1 ari ty  laws, 
which  do  not  need  this  information. 


The  following  particular  form  of  the  periodic  free-stream 
velocity  is  now  assumed: 


U„(t)  =  ,J0  (1+G  sin  wt)  (4.5) 

where  both  the  time-mean  value,  UQ>  and  the  relative  amplitude  of 
oscillation,  e,  are  constants.  Note  that,  from  Eq.  (4.5),  Ue  =  u0, 
i.e.,  the  time-mean  flow  is  a  "flat  plate"  type  of  boundary  layer,  in 
the  sense  that  it  is  not  subjected  to  any  effective  pressure 
gradient.  In  fact,  the  time-mean  quantities  are  related  through  the 
following  equations: 


9U 

3X 


IT  M  +  w  M  _  82U 

u  3X  +  V  3y  -  v  -^2  +  3y 


-  3U. 


3__U  +  ll  _  '  P 


3x 


3U  V 
P  P 


3y 


(4.6) 

(4.7) 


and  boundary  conditions 


U  =  0,  V  =  0 

at  y  =  0 

u  ♦  u 

0 

as  y  +  °° 

(4.8) 

^  1 

II 

o 

at  y  =  0 
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resulting  from  time-averaging  Eqs.  (4.1)  to  (4.4).  The  last  two 

terms  of  Eq.  (4.7)  are  the  Reynolds  stresses  due  to  the  oscillatory 

2 

motion.  They  are  of  the  order  e  ,  and  can  be  neglected  provided  the 
amplitude  of  oscillation  of  the  freestream  is  sufficiently  small, 
i.e.,  e  <<  1.  Assuming  this  to  be  the  case,  the  above  system  reduces 
to  the  standard  boundary-1 ayer  equations  for  a  flat  plate.  However, 
differences  with  respect  to  the  correspondi ng  steady  motion  (i.e., 
the  one  produced  by  a  free-stream  velocity  Ue  =  u0)  could,  in 
principle,  arise  if  the  time-mean  turbulent  shear  stress,  x,  is 
affected  by  unsteadiness.  It  has  been  speculated  that  this  could 
happen  if  the  forcing  frequency  f  =  w/2tt  is  high  enough  to  interact 
with  the  most  energetic  turbulent  frequencies  [Ramaprian  and  Tu 
(1982)].  A  measure  of  this  frequency  range  is  provided  by  the 
bursting  frequency,  fp,  which  in  a  boundary  layer  can  be  correlated 
with  global  quantities  through  the  approximate  relation  [Rao  et  al. 
(1971)] 


U 


Ramaprian  and  Tu  (1982)  give  some  experimental  evidence.  Past  and 
present  experiments  have  indicated  that,  in  any  case,  the  resulting 
variation  in  the  time-mean  velocity  profile  relative  to  the 
correspondi  ng  steady  profile  is  small  for  large  Reynolds  numbers.  It 
can,  therefore,  be  neglected  for  purposes  of  the  perturbation 
analysis  presented  here. 
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The  assumption  e  <<  1  allows  a  linear  analysis  in  the  amplitude 
to  be  performed.  This  means  that  only  motions  at  the  forcing 
frequency  are  considered  to  be  significant,  i.e.,  higher  harmonic 
oscillations  can  be  neglected.  The  advantage  of  the  complex  notation 
can,  then,  be  utilized,  defining 


U  =  U  +  eU1  e 


iut 


V  =  V  +  E  V, 


i  cut 


T  =  T  +  £  e 


iut 


and  restating  Eq.  (4.5)  as 


ue  =  V1  +  e  eiwt) 


(4.10) 


(4.11) 

(4.12) 


(4.13) 


Note  that  U^,  and  are  all,  in  general,  complex  quantities. 
Introducing  Eqs.  (4.10)  to  (4.13)  into  Eqs.  (4.1)  to  (4.4),  and  using 
Eqs.  (4.6)  to  (4.8),  the  equations  and  boundary  conditions  for  the 
oscillatory  components  are  obtained  as  [Telionis  (1981)] 


3U1  9V. 

sT  *  "sy  '  0 


(4.14) 


-  3U1  3U  _  9Ui  77 

f»(uru0)  *  U1  «  tV7r+  V1  w 


=  v 


92U1  9tj_ 


ay 


2  ay 


(4.15) 


t.'f. 


u,  -*■  u 
1  0 


as  y  > 


at  y  =  0 


(4.16) 


Before  undertaking  the  analysis  of  Eqs.  (4.14)  to  (4.16),  the 
asymptotic  behavior  of  the  time-mean  flow  for  large  Reynolds  number 
is  briefly  discussed,  in  the  next  section. 


4.3.  Asymptotic  Characteri sti cs  of  the 
Time-Mean  Flow  at  Large  Reynolds  Numbers 

As  pointed  out  in  the  previous  section,  the  time-mean  velocity 
is  practically  the  same  as  for  the  corresponding  steady  flow,  to  the 
first  order  in  the  amplitude  of  oscillation.  Tennekes  and  Lumley 
(1980)  presented  an  asymptotic  analysis  of  the  steady  flow  for  very 
large  Reynolds  numbers,  without  using  any  closure  hypothesis  for  the 
turbulent  shear  stress.  Their  results,  then,  can  be  assumed  to  apply 
also  to  the  time-mean  unsteady  flow.  Their  main  conclusions,  with 
reference  to  the  present  context,  can  be  summarized  as  follows: 

(i)  There  exist  two  distinctly  different  length  scales  for  the 
motion,  which  determine  a  double-1 ayei  structure:  an  outer  layer 

characteri zed  by  the  (local)  boundary  layer  thickness,  5,  and  an 
inner  layer  with  length  scale  v/u#.  The  quantities  S(x)  and  u+(x) 
will  be  assumed  to  correspond  to  the  steady  flow  at  the  time-mean 
velocity.  41  ternati vely ,  they  could  be  chosen  as  the  time-mean 

values  in  unsteady  flow  which  cannot  differ  much  from  the  previous 


rti 


ones  due  to  the  small  amplitude  of  oscillation.  Furthermore,  as  the 
boundary  layer  thickness,  6,  is  an  ill-defined  quantity,  a  thickness, 
A,  defined  by 


will  be  used  instead.  Mote  that  the  ratio  of  the  two  length  scales 
defines  the  following  Reynolds  number: 


(4.18) 


With  the  present  definition  of  a,  Eq.  (4.17),  it  follows  that  R* 
=  Re  =  U  0/v.  However,  it  is  better  to  work  with  R*  in  order  to  be 

U  0 

able  to  generalize  the  results  to  fully  developed  flow  in  pipes  and 
two-dimensional  channels  (see  section  4.6). 

(ii)  The  velocity  scale  for  both  the  layers  is  u*.  Then,  the 

2 

time  scales  are  tQ  =  a/u*  and  t^  =  v/u*  for  the  outer  and  inner 
layers,  respectively. 

(iii)  In  the  limit  R*  +  »,  the  following  relations  hold: 

U-U 

“u —  =  fr('^)  for  the  outer  layer  (4.19) 

7j  yu* 

—  =  f  (— — )  for  the  inner  layer  (4.20) 

where  the  functions  F  and  f  can  only  be  obtained  from  experiments. 
This  indeterminacy  arises  because  the  system  of  equations  is  not 
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closed.  Eqs.  (4.19)  and  (4.20)  are  the  well-known  "vel oci ty -defect 
law"  and  "law  of  the  wall". 

(iv)  For  R*  +  °°,  there  exists  a  region  where  the  outer  and 
inner  layers  overlap.  In  this  "inertial  sublayer"  Eqs.  (4.19)  and 
(4.20)  are  both  valid,  leading  to  the  following  results 

TT  1  yu* 

U-tl  . 

,  I  (i)  .  D  (4.22) 

where  k  (von  Karman  constant),  B  and  !)  are  constant  and  independent 
of  the  Reynolds  number  R*.  Equation  (4.21)  or  (4.22)  is  the  well- 
known  "logarithmic  law". 

(v)  Subtracting  Eq.  (4.22)  from  Eq .  (4.21)  yields 

1  U  1 

^  ^  In  R*  +  B-0  (4.23) 


which  is  called  the  "logarithmic  friction  law",  as  it  allows  the 
determination  of  the  shear  velocity  (and,  then,  the  wall  shear 
stress)  if  the  constants  k,  B  and  D  and  the  boundary  layer 
thickness  a  are  known. 

(vi  )  If  L  is  the  length  scale  in  the  longitudinal  direction, 

then  a/L  0  as  R+  >  ».  Equation  (4.23)  also  shows  that 

E  +  Q  as  R^  -*■  4°.  Howe /or,  the  ratio  _  u  L,,j  A  remains  finite  in 

*  '  '  0" 

the  limit.  The  length  scale  t  can  he  correlated  to  local  quantities 
by  performing  an  order-of-ingo’ t  ide  mnlysis  in  the  moneotum  integral 
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equation,  which  for  a  steady  flow  at  zero  pressure  gradient  is  simply 
[White  (1974)] 


d9  p2 
dx  =  L 


(4.24) 


Then,  L  can  be  defined  as 


(4.25) 


Notice  that,  from  Eqs.  (4.17)  and  (4.25),  it  follows  that  y  =  1,  and 
hence  the  definition  of  L  is  consistent  with  the  asymptotic 
requirement  that  y  should  remain  finite. 

The  foregoing  results  will  be  very  usef  in  the  analysis  of  the 
oscillatory  motion,  discussed  in  the  following  section. 


4.4.  Asymptotic  Characteristics  of  the 
Oscillatory  Motion 

The  imposed  oscillation  introduces  a  third  time  scale,  into 
the  problem.  Its  relation  to  the  outer  and  inner  time  scales  can  be 
described  by  the  following  two  frequency  parameters 


~  o  wA  wL 

"  =  -1  =  u7  =  TT 

w  *  o 


(4.26) 
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In  addition,  a  third  length  scale,  namely  u*/u>,  appears,  being 
associated  with  the  imposed  oscillation.  This  defines  an  "unsteady 
layer"  whose  thickness  bears  to  the  outer  and  inner  thicknesses  the 
same  relation  as  the  corresponding  time  scales,  and  is  thus  described 
by  the  same  parameters,  i.e.,  Z  and  a,  defined  above. 

Various  flow  regimes  can  be  distinguished,  according  to  the 
values  of  the  foregoing  parameters.  For  Z  <  1,  one  has  that 

t  and  hence,  significant  unsteady  effects  must  be  expected  in 
~  0 

the  outer  layer.  In  fact,  the  unsteady  layer  is  merged  with  the 

outer  layer.  On  the  contrary,  a  =  R*/w  >  R*  >>  1,  i.e.,  m  >>  t. 
indicates  that  the  inner  layer  responds  instantly  to  the  excitation, 
i.e.,  in  a  quasi-steady  manner.  This  will  be  called  the  "low 
frequency  regime".  Note  that,  in  the  particular  case,  to  <<  1,  the 
outer  layer  will  also  show  an  instantaneous  response,  i.».,  the 
entire  boundary  layer  will  behave  in  a  quasi -steady  manner, 

The  other  extreme  situation  arises  when  ^  >  R*  i.e.  a  <  1,  in 
which  case  <  t^.  In  this  case,  the  unsteady  layer  will  be 

actually  submerged  in  the  inner  layer.  The  effects  of  unsteadiness 
are,  therefore,  significant  in  the  inner  layer.  On  the  other  hand, 
since  <<  t  ,  there  is  not  enough  time  for  the  turbulent  flow  in 
the  outer  layer  to  respond  to  the  imposed  oscillation  and  hence  this 
layer  will  move  as  a  rigid  body,  i.e.,  it  will  perform  inviscid  or 
slug-flow  osci  1 1  ations.  In  fact,  even  the  outer  part  of  the  inner 
layer  may  exhibit  the  slug-flow  behavior.  This  situation  will  be 
named  "very  high  frequency  regime". 


142 


The  situation  1  <<  u  <<  R*  which  lies  between  the  two  1  i  mi  tiny 

cases,  described  above,  needs  careful  study.  Since,  under  this 

condition,  the  inner  and  outer  time  scales  are  such  that  t^ 

<<  w”*  «  tQ,  one  should  expect  that  there  will  be  quasi -steady  flow 

in  the  inner  layer,  while  a  slug-flow  is  present  in  the  outer 

layer.  Furthermore,  the  fact  that  v/u*  <<  u^/u  <<  a  in  this  flow, 

suggest  that  the  unsteady  layer  lies  between  the  outer  and  inner 

layers.  Only  in  this  unsteady  layer,  then,  the  unsteady  viscous 

effects  are  significant.  In  fact,  it  has  been  found  to  be  useful  in 

the  analysis  to  subdivide  this  regime  further  into  two  distinct 

regimes:  an  "intermediate  frequency  regime",  for  which 

~  1/2  1/2 

1  <<  a)  <  R*  (or,  equivalently,  R*  <  a  <<  R*),  and  a  "high  fre- 

1/2  ~ 

quency  regime",  with  Rj  <<  w  <<  R*(or  ,  equivalently,  1  <<  a 

«  ^y2).  Notice  that  u>  =  corresponds  to  the  case  u>_1  = 

(toti)1^2*  i*e*  the  geometric  mean  of  the  outer  and  inner  time 
scales,  for  which  Z  =  a. 

Each  one  of  these  four  frequency  regimes  is  analyzed  in  detail 
in  the  following  sub-sections. 

4.4.1.  Low  Frequency  Regime 

This  regime,  as  explained  above,  is  characterized  by  the 
conditions  w  <  1,  a  >  R*.  The  outer  layer,  where  unsteady  effects 
are  expected  to  be  important,  is  analyzed  first.  For  this  purpose, 
the  following  dimensionless  variables  are  introduced: 


w 


The  quantity  F  defined  in  Eq.  (4.28)  is  the  same  as  in  Eq.  (4.19), 

and  hence  F  =  F(n).  The  definition  of  the  variable  ,  in  turn,  stems 

from  the  similarity  between  the  equations  and  boundary  conditions 

for  U  and  U^.  In  the  other  hand,  the  nondimensional  ization  of  the 

transverse  velocity  components  in  Eq.  (4.28)  is  suggested  by  the 

respective  continuity  equations.  Finally,  the  scale  u*2  for  the 

oscillatory  component  of  the  turbulent  shear  stress,  ip  arises  from 

the  following  quasi -steady  argument.  As  already  mentioned  in  Chapter 

3,  it  is  well  known,  from  experiments  in  steady  flow, 

2 

that  x/Ue  depends  only  weakly  on  Ue  (i  .e. ,  on  the  Reynolds 

number).  Considering  it  to  remain  constant,  to  a  first 

approximation,  when  Ue  is  varied  by  All  ,  the  variation  At  in  t  is 

related  to  aU  by 
e 


(4.29) 


Equation  (4.29),  together  with  Eq.  (4.12)  and  (4.13),  show 
that  ^  =  0(ti*)  in  quasi -steady  flow.  As  shown  in  chapter  3, 

experiments  however  suggest  that  u*2  i  n  fact,  the  appropriate 


scale  also  for  any  finite  frequency. 


Introducing  the  definitions  (4.28)  into  Fqs.  (4.14)  and  (4 
one  obtains 


where 
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In  the  limit  +  °°,  Eq.  (4.31)  reduces,  to  the  lowest  order,  to 


3°1  9flfl  dF 
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The  outer  boundary  condition  for  Ui  in  Eq.  (4.16),  can  be  expressed 


Uj  +  0 


as  n  ® 


^  H  .  J  3  j 


According  to  Eq.  (4.34),  viscous  effects  on  the  oscillatory  motion 
are  negligible  in  the  outer  layer,  to  this  order  of  approximation. 
One  the  other  hand,  from  both  Eqs.  (4.30)  and  (4.34)  it  is  observed 
that,  for  similarity  solutions  to  exist,  i.e.,  solutions  independent 
of  the  streamwise  coordinate  5, it  is  necessary  that  the 
coefficients  and  ^  be  constant.  Not  surprisingly,  these  are  the 
same  conditions  as  those  necessary  to  achieve  self-preservation  in 
the  time-mean  flow  [Yajnik  (1970)].  Now,  from  Eqs.  (4.17),  (4.24), 
and  (4.25),  it  can  easily  be  shown  that  p^  =  1  -  p_.  Hence,  the  only 
necessary  condition  for  similarity  is 


90  2  du* 

3~  dx“  =  const* 
u* 


(4.36) 


Now,  for  very  high  Reynolds  numbers  one  has  that  u^/U^  =  const.  This 
shows  that  Eq.  (4.36)  is  satisfied,  with  the  constant  on  the  right 
hand  side  equal  to  zero,  and  that  similarity  is,  in  fact,  possible  in 
this  case.  Then,  a  "velocity -defect  law"  for  the  oscillatory  motion 
can  be  expected  to  hold  in  the  form 


U.-U 

»i 2  -uT  ■ F.  => 


(4.37) 


where  -j  s  a  complex  function  that  must  be  determined  from 

experiments.  Note  that  Eq.  (4.37)  is  riot  a  sel  f  -preservi  ng  solution 
for  any  physical  (constant  <o)  oscillatory  flow,  because  still 


depends  on  x  (or  r ).  Self-si  ri  ’  in  ty  in  a  given  f!  ow  can  be  achieved 


only  in  the  quasi -steady  limit,  i.e.,  for  ]>  <<  1,  for  which  the 
effects  of  the  local  acceleration  and  the  pressure  gradient  [lumped 
together  in  the  first  term  of  Eq.  (4.34)]  are  negligible.  In  this 
case,  it  can  be  shown  that  F^(^)  is  practically  coincident 
with  F(n).  However,  Eq.  (4.37)  should  still  be  useful  in  providing 
the  link  between  oscillatory  motions  at  different  frequencies. 

Next,  the  inner  layer  is  analyzed.  The  nondimeosi onal i zati on  Is 
carried  out  using  similar  criteria  as  for  the  outer  layer,  with  the 
resul ts 


rl 


(4. 38) 


Note  that  the  function  f  in  Fq.  (4.38)  is  the  same  as  in  Eq.  (4. 2D); 
then,  f  =  f(y+).  Introducing  the  definitions  (4.38)  into  Eqs.  (4,14) 
and  (4.15)  yields 
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When  R.  «,  the  zero  order  approximation  of  Eq.  (4.40)  is  simply 


A:  a,; 

+7  +  +  s  0 

sy  3y 


(4.41) 


i.e.,  a  balance  between  viscous  and  turbulent  diffusion.  The  inner 
boundary  conditions  of  Eq.  (4.16),  are  now  rewritten  as 


=  0,  V, 


(4.42) 


Note  that  (as  expected)  Eq.  (4.41)  will  yield  quasi -steady  solutions, 
since  it  does  not  depend  on  w.  Furthermore,  sel  f  -simi  1  ari  ty  can  be 
expected  to  hold  in  the  form  of  a  "law  of  the  wall",  i.e. 


fi 


(4.43) 


with  the  function  fj  to  be  determined. 

Now,  as  in  the  case  of  the  time-mean  flow,  an  asymptotic 
matching  procedure  can  be  applied  to  the  inner  and  outer  solutions 
for  the  oscillatory  flow.  The  procedure  proposed  by  Van  ilyte  (1978) 
will  be  followed  here.  Symbol i cal ly ,  it  states  that 


E 


E 


In  words,  Eq.  (4.44)  expresses  that  the  m-term  outer  expansion  of  the 
n-tenn  inner  expansion  of  the  quantity  between  brackets  equals  the  n- 
tenn  inner  expansion  of  its  m-term  outer  expansion.  In  the  present 
case,  the  expansions  are  in  terns  of  the  large  parameter  R*,  and  the 
quantity  chosen  to  be  matched  is  u|  =  U,/u*.  Now,  as  in  the  inner 
layer  u|  is  assumed  to  satisfy  Eq.  (4.43),  it  does  not  have  any 
explicit  dependence  on  R*.  Hence  the  inner  expansion  of  u|  consists 
of  only  one  term,  namely  fj(y  +  ).  On  the  other  hand,  from  Eqs. 
(4.23),  (4,37)  and  (4.38),  one  gets  the  following  outer  expansion 


(n;  R*)  =  £•  In  R*  +  [B-D  +  F1  (n  ?«)]  (4.45) 


i.e.,  a  two  term-expansion,  with  the  leading  term  growing  very  large 
for  R*  ->•  «>.  The  matching  procedure  is  started  with  m  =  n  =  1.  Then, 


I!°1(ul)  =  Mf  ln  R*>  =  F  ln  R* 


(4.46) 


01I1(Uj)  =  01[f1(R*n)] 


(4.47) 


From  the  matching  principle,  Eq.  (4.44),  and  from  Eqs.  (4.46)  and 
(4.47),  one  obtains 


0,  [f  i  (R*n)]  =  T7  ^  R, 


(4.48) 
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which  is  satisfied  only  if 
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f  1  (y+ )  ^  In  y+  +  B1  as  y+  +  ®  (4.49) 

where  is  an  unknown  constant.  Eq.  (4.49)  shows  that  the  law  of 

the  wall  merges  with  a  logarithmic  law  when  y+  is  large.  The  next 
step  in  the  matching  is  to  take  m  =  2  and  n  =  1.  Then, 

Ii°2(ui)  =  ln  R* +  B-D  +  Fi(iv  ^)] 

,  + 

=  7  In  R*  +  I}  [Fj  (^;  S)]  (4.50) 

°?.h  (,J1>  =  02[fi(R*n)] 

=  IT  >"  (R*n)  +  Bj  (4.51) 

Equation  (4.49)  has  been  used  in  Eq.  (4.51).  From  Eqs.  (4.44), 
(4.50)  and  (4.51)  one  obtains 

i1cf1(r7;  ^)3  =  r  ln  n  +  Bi  (4.5?) 


which  means  that 
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where  0^  is  a  new  constant.  From  Eq.  (4.53)  it  is  seen  that  the 
velocity  defect  law  also  merges  with  a  logarithmic  law  when  n  0. 

Interestingly,  Eqs.  (4.49)  and  (4.53)  show  that  the  slope  of  the 
logarithmic  curve,  1/k,  is  the  same  as  for  the  time-mean  flow. 

4.4.2.  Very  High  Frequency  Regime 

The  conditions  satisfied  by  the  frequency  parameters  are  now 
to  >  R*,  a  <  1.  The  lowest  order  approximation  of  Eq.  (4.31)  is, 
then,  simply 


(4.54) 


i.e.,  a  balance  between  local  acceleration  and  pressure  gradient. 
Equation  (4.54)  shows  that  the  outer  layer  performs  rigid  body 
oscillations,  as  expected. 

On  the  other  hand,  when  the  zero  order  approximation  is  sought 
in  the  inner  momentum  equation,  Eq.  (4.40),  a  contradiction  is  found, 
since  the  leading  term  is  a  constant.  This  indicates  that  there 
exists  a  scaling  problem.  In  fact,  Eq.  (4.54)  suggests  that  changes 
in  the  velocity  of  the  order  of  U0  must  occur  within  the  inner  layer 
in  order  to  match  the  outer  solution.  Hence,  the  following  alterna¬ 
tive  nondimensional i zation  is  proposed: 


(4.56) 
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after  which,  the  zero-order  approximation  of  Eq.  (4.40)  becomes 

2  - . + + 

,,  3  U, 

iul  TT  (4‘57) 

ay 

i  .e.  a  balance  between  local  acceleration,  pressure  gradient  and 
viscous  diffusion.  In  turn,  the  inner  boundary  condition  for  the 
velocity  is,  from  Eqs.  (4.42)  and  (4.55), 

ll|+  =  -  1  at  y  =  0  (4.58) 

while  the  matching  condition  with  the  outer  layer  simply  becomes, 
from  Eg.  (4.54) 

*  0  as  yT  +  -  (4.59) 

The  solution  to  Eqs.  (4.57)  to  (4.59)  is  the  well-known  Stokes  or 
shear-wave  solution,  namely 

ut*  =  -  exp  f-  (4.60) 

1  '  /7a 

which  is  obviously  self-similar  in  the  inner  coordinates.  Hence,  for 
these  very  high  frequencies,  the  oscillatory  motion  behaves  as  in 
laminar  flow  [Lighthi 1 1  (1954)].  Note,  from  Eq.  (4.60),  that  the 

main  change  in  the  oscillatory  velocity  occurs  in  a  layer  whose  width 
is  of  the  order  of  xs,  where 


is  the  Stokes  thickness.  For  a.  ~  1,  the  Stokes  layer  thickness  is  of 
the  same  order  as  the  inner  layer  thickness.  However,  for 
a  <<  1  (w  >>  R*),  6S  becomes  much  smaller  than  the  thickness  of  the 
inner  layer,  which  will  then  perform  slug-like  oscillations,  just 
like  the  outer  layer. 


4.4.3.  Intermediate  Frequency  Regime 

The  conditions  that  characterize  this  frequency  regime  are 
~  1/2 

1  <<  w  <  R*  <  a  <<  R*.  Hence,  the  zero  order  approximation  in  both 
R*  and  u  of  the  outer  momentum  equation,  Eq.  (4.31),  is  again  Eq. 
(4.54),  i.e.  there  exists  an  outer  slug  flow,  as  in  the  very  high 
frequency  regime. 

1/2 

In  turn,  taking  into  account  that  aE  >  (R*  /In  R*) 

+  °°,  as  R+  +  od,  Eq.  (4.40)  for  the  inner  layer  reduces  to  Eq. 
(4.41).  Hence,  there  is  a  quasi-steady  motion  in  the  inner  layer, 
and  a  1 aw-of -the-wal  1  like  Eq.  (4.43)  can  be  expected  to  hold.  Note 
that,  this  time,  nothing  suggests  that  the  scaling  for  the 
oscillatory  velocity  is  wrong:  neither  the  momentum  equation  nor  the 
matching  with  the  outer  layer,  as  this  is  performed  through  the 
unsteady  layer,  analyzed  next. 

The  dimensionless  variables  for  the  unsteady  layer  are  defined 
as  follows: 


Note  that  F  now  depends  both  on  y  and  The  scaling  for  Uj  seems  to 
be  appropriate  to  bridge  the  gap  between  the  inner  and  outer  flows. 
Introducing  the  definitions  (4.62)  into  Eqs.  (4.14)  and  (4.15)  one 
obtains 
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To  zero  order  in  to  \  Eq.  (4.64)  reduces  tr 


i'J, 


3t, 


ay 


(4.66) 


i.e.,  a  balance  between  local  accelerat ion,  pressure  gradient  and 
turbulent  diffusion.  E  |.  (4.65)  ad:iifs  s'*1  ‘  •  *.  ini  1  ar  solutions  in  the 


fora  of  a  velocity  def 


t  law,  i.e. 
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(4.56) 


with  F2  to  be  determined.  Matching  with  the  outer  layer  simply 
requi res 


Fo(y)  +0  as  y  +  (4.67) 

On  the  other  hand,  matching  with  the  inner  layer,  using  Van  Dyke's 

principle,  is  more  cumbersome  than  in  the  low-frequency  regime,  as 

the  three  parameters  w,  a  and  R*  are  now  assumed  to  grow  very 

large.  For  1  «  u  «  R*  ,  however,  u  may  be  considered  fixed 

when  R*  >  and  the  matching  principle  indicates,  as  before,  the 

1/2 

existence  of  a  logarithmic  region.  As  u>  -►  R*  ,  this  region  would 

tend  to  disappear. 

4.4.4.  High  Frequency  Regime 

This  regime  is  characteri zed  by  the  conditions  1  <<  a  <<  R*^ 
<<  a!  <<  R*.  The  outer  layer  motion,  hence,  continues  to  be  a  slug 
flow  oscillation,  as  Eq.  (4.31)  reduces  to  Eq.  (4.54). 

'  In  the  inner  layer,  Eq.  (4.40)  leads  again  to  a  contradiction, 
as  aE  +  0  when  R*  +  «  (a  is  now  fixed  when  R+  *  ®).  Using  the 

alternate  scaling  proposed  for  the  very  high  frequency  regime,  Eqs. 

(4.55)  and  (4.56),  the  zero-order  approximation  of  Eq.  (4.40)  is  then 
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(4.68) 


Eq.  (4.68)  shows  that  the  only  significant  mechanism  is  now  viscous 
diffusion,  which  must  be  counteracting  the  difference  in  shear  stress 
at  the  two  ends  of  the  inner  layer,  i.e.,  it  is  undergoing  a  Couette 
motion.  From  Eq.  (4.68)  an  I  the  inner  boundary  condition,  Eq. 
(4.58),  the  following  solution  is  obtained 

' 1  j !  -  a /  *  -  1  (4.69) 

where  a  is  a  constant  that  depends  on  the  conditions  in  the  unsteady 
1  aver. 

The  analysis  of  the  unsteady  layer  does  not  differ  at  all  from 
the  one  for  the  intermediate  frequency -regime  case.  Thus,  Eqs. 
(4.65)  to  (4.67)  are  still  expected  to  nold.  Mahhing  with  the  inner 
layer,  in  turn,  requires  that  F?  become  a  linear  function  of  its 
argument  when  this  becomes  very  small.  Using  Eg.  (4.69),  this  can  be 
written  as 

F2(y)  -  as  y  -►  0  (4.70) 

4.4.5.  Summary  of  Results 

The  main  results  of  the  foregoing  analysis  are  now  summarized. 
The  original  concept  introduced  is  that  of  a  triple-layered  structure 
imposed  by  unsteadiness,  as  schematized  in  Fig.  4.1(b).  This  is 
superimposed  oi  the  original  two-1  oyer  ruotw'e,  already  present  in 
steady  flow,  shown  in  Fig.  4.1(a).  The  different  flow  regimes  are 


defined  according  to  the  way  the  unsteady  structure  overlaps  the 
steady  structure.  This  is  schenatical ly  represented  in  Fig.  4.2. 
Different  similarity  laws  hold  in  the  different  layers,  for  each 
frequency  regime.  These  laws  were  derived  in  terms  of  the  complex 
quantity  U},  defined  by  Eq.  (4.10).  They  are  restated  below  in  a 
form,  that  is  more  convenient  to  use  in  practice. 

For  a  free-stream  velocity  given  by  Eq.  (4.5),  the  longitudinal 
velocity  component  given  by  Eq.  (4.10)  is  written  as 


U  =  U  +  e  ('J,  j  sin  wt  -  IJ^  COSiot)  (4.71) 

where  U-q  and  are,  respectively,  the  in-phase  and  out-of-phase 
components  of  the  oscillatory  motion.  To  obtain  the  results  for 
and  IJ^2  from  those  for  Uy,  it  is  only  necessary  to  replace  Uj  by  ('J^ 
-  UI12)  in  the  various  relations  derived  earlier.  The  following  are 
the  results  for  and  2- 


(i)  Low  frequency  regime:  Z  <  1,  a  >  R* 

(a)  Outer  layer  (coincident  with  unsteady  layer) 
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(b)  Unsteady  layer: 


Vs 
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(4.37) 


where  a^  and  a?  are  unknown  constants, 


(iv)  Very  high  frequency  regime:  Z  >  R*,  «  <  1 


(a)  Outer  layer;  same  results  as  for  the  intermediate 
frequency  regime. 


(h)  Inner  layer  (coincident  with  unsteady  layer) 
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(4.89) 


The  foregoing  results  are  asymptotic  limits  for  -»■  ®.  Hence,  they 
should  be  referred  to  as  "weax"  similarity  laws  [YajniK  (1970)],  in 
the  sense  that  they  hold  only  to  the  lowest  order  asymptotic 
analysis.  Deviations  from  these  asymptotic  laws  should  be  expected 
in  practice,  due  to  the  finite  value  of  R*. 


4.3.  Generalization  to  Nonzero  Time-Mean 


Pressure  Gradient  Round ary  Layer 

The  concents  developed  so  Tar  for  zero  time-mean  pressure 


gradient  flows  admit  a  relatively 
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flow.  Not  onl  /  does  such  an  •>  •  tens ;  w>  : 
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application,  but  also  enables  one  to  verify  some  of  the  general 


conclusions  using  the  available  experimental  i nformation.  In  this 
section,  the  general  case  of  nonzero  time-mean  pressure  gradient 
boundary  layers  is  examined.  Fully  developed  channel  and  pipe  flows 
are  analyzed  in  the  next  section. 

Consider  a  general  harmonic  free-stream  oscillation  of  the  form 


Ug(x,t)  =  'j^(x)  [1  +  e (x )  sin  ^t]  (4.90) 


with  c  «  1.  To  the  first  order  in  the  relative  amplitude,  z,  the 
driving  pressure  gradient  is,  from  Eqs.  (4.3)  and  (4.90), 
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(4.91) 


The  two  extra  contributions,  namely  the  first  and  the  last  terms  on 
the  right  hand  side  of  Eq.  (4.91),  will  produce  modifications  in  the 
time-mean  flow  and  the  oscillatory  motion,  respectively,  relative  to 
the  case  of  zero  time-mean  pressure  gradient.  Note  that,  rigorously 
speaking,  flows  with  =  const.,  c  =  c(x)  would  belong  to  the 
category  of  zero  pressure-gradient  flows.  However,  this  name  will  be 
reserved  for  the  situation  in  which  both  l)  and  .  are  constant. 


16 1 


4.5.1.  Time -Moan  Flow 

The  discussion  on  the  features  of  the  1 1 me -mean  flow  of  section 
4.3  can  be  generalised  to  include  the  present  case  in  the  following 
manner  [Tennekes  and  Lumley  ( ."')) 


(i)  The  normalised  houwdury  layer  thickness,  used  as  the 
length  scale  for  the  outer  layer  is  redefined  as 
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(4.92) 


which  is  a  more  general  expression  than  Kg.  (4.17). 


(ii)  In  addition  to  the  out'-r  and  nner  time  scales,  the 
longi tudina?  variation  if  the  free -stream  velocity 
introduces  its  own  t i  scale,  namely  j -rJU^/d k (  *^.  Its 
ratios  to  the  preceding  time  scales  define,  respect i vely , 
the  following  two  parameters : 
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infinitely  large  compare.-*,  *.o  the  scale  of  t"e  inner  layer 
at  this  order  of  approximation. 

(iii)  When  R+  >  -*>,  the  vel  or  i  ty  -defer*  law.  Eg.  (4.19), 
generalizes  to 

v  * F  if  o'  !,-55) 

For  Eq.  (4.95)  to  hold,  it  is  required  that  dy./dx  =  0. 
This  condition  defines  the  so  called  equilibrium  boundary 
layers,  for  which  Eq.  (4.95)  gives  a  sel f -preservi ng 
velocity  profile  (in  the  outer  layer).  On  the  other  hand, 
the  law  of  the  wall,  Eq.  (4.20),  is  still  valid. 

(iv)  The  logarithmic  laws,  Eqs.  (4.21)  and  (4.22),  are 
satisfied  as  long  as  remains  finite  [see  (vi )  below]. 
However,  the  extent  of  the  inertial  sublayer  and  the  value 
of  0  in  Eq.  (4.22)  depend  on 

(v)  A  logarithmic  friction  law  like  that  of  Eq.  (4.23)  holds 
for  finite 


(vi)  The  momentum  integral  e.-g.i,i  t i O"  for  steady  flow  is  [White 
(1974)] 
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Equation  (4.9b)  shows  that,  in  addition  to  the 
1  oriqi  tudi  na  1  length  scale  I  defined  b.  :  j.  (4.25),  a 
second  length  suM,  i  M  /vd"  /dx),  is  Moused  by  the 
spatial  variation  of  t!u;  free  stream.  It  will  to  assumed 
that  i  >  1,  i  ,e.  rapidly  accelerat inq  or  decelerating 
flows  are  excluded.  Hi  on.  I'd  /dx  <  11  /!.,  as  assumed 

•:o  -  »a) 

above,  and  <  1. 


4.5.2.  Oscillatory  Motion 

Regarding  the  oscillatory  iiiotiou,  the  discussion  presented  at 
the  beginning  of  section  4.4  still  applies,  if  in  Eq.  (4.2(1)  U0  is 
changed  to  M  .  Retaining  tb.-  definitions  (4.  Id'  to  (4.12)  and 

CO 

updating  Eg.  (4.13)  to 

i)  ,  !  Ms  o'  r'l  (1.97) 

Equation  (4.15)  can  now  !v>  geni-ra  1  i  :>i  \  as 


The  outer  seal  ini] 


■it  F/|.  { 4 .  f?  >< )  is  still  -appropriate  if  Jt>  is 
replaced  by  ij  .  Then,  vi  e  ••:..)  res  i  on  an  3 ’onons  t;  Eq.  (4.31)  results, 

CO 

except  for  the  following  a 'ldi  tinea  I  term  on  its  left  hand  $  id-* 

-  (•■•1  "  19)  (’J j_  +  f) 

where 

t  d_q 

r2  r:  *  "e  dx 

is  a  new  "equilibrium"  parameter.  In  the  low  frequency  regime,  then, 
in  order  for  a  similarity  law  like  Eq.  (4.37)  to  hold  (with  U0 
replaced  by  U  ) ,  it  is  necessary  that,  in  addition  to  satisfying  Eq. 

00 

(4.36),  both  and  be  constant.  For  the  remaining  flow  regimes, 
in  turn,  u  >>  1  and,  hence,  the  flow  will  be  slug-like  in  the  outer 
layer,  as  for  zero-pressure-gradient  flows. 

The  analysis  for  the  inner  layer  remains  exactly  the  same  as 
before  because,  as  explained  above,  r+  [Eq.  (4.94)]  vanishes  at  the 
lowest  order  of  approximation.  Finally,  the  equation  for  the 
unsteady  layer,  in  the  intermediate  and  high  frequency  regimes,  is 
still  analogous  to  Eq.  (4.64),  except  for  the  following  additional 
term  on  its  left  hand  side 


(4.99) 

(4. LOO) 
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where  the  same  scaling  as  before  is  us-'d,  name ly  tq.  (4.62),  with  U0 
replaced  by  II  .  q  v,.  j  for  boih  regimes,  Fq.  (4.65)  is  again 

valid  to  Imie  lowest  order  of  appro ■  in>  1 1  i on  and.  1  m  a  similarity 


law  like  that  in  F.q.  (4.66)  can  im  expect. *d  to  hold  (with  !l0  replaced 
by  U  ). 

oo 

Though  the  above  discussion  was  limited  to  equilibrium  boundary 
layers,  only  a  few  of  the  result*:  ao  a:.tua!|v  t  rir  ted  by  this 
condition.  First,  for  nonegui  I  i  bri  urn  bom  tar/  layer..,  the  velocity- 
defect  law,  Fq.  (4.96',  is  not  valid,  as  the  profile  will  also  depend 
on  its  "history",  i.e.,  its  upstream  evolution.  However,  a 
logarithmic  region  still  exists  (for  finite  y1 )  though  its  extent 
(and  the  value  of  II)  will  depend  on  the  history.  Then,  the  only 
difference  in  the  results  for  the  oscillatory  motion  is  that  no 
similarity  law,  like  that  of  Fq.  (4,37),  wi  1  i  hold  in  the  outer  layer 
for  the  low  frequency  regime. 


4.6.  Extension  totally  Dove  lope  !  Channel 
and  Pipe  Flows 

In  fully  developed  flow  all  the  properties  are  invariant  in  the 
longitudinal  direction.  In  particular,  the  shear  velocity,  u*,  is  a 
constant  (in  the  corresponding  steady  flow).  Also,  since  the 
convection  terns  vanish  (as  d  m  .  the  cross ->trea..i  velocity  <V>)  the 
problem  becomes  linear  and  the  results  are  valid  for  any  amn1 i tude  of 
oscillation,  not  necessan  ] nil,  gro/ided  l  instantaneous 
Reynolds  number  re.su  '  n  s  !  irge  f-.i'n  ,  /i:,o  •  ne  ,  j  le. 


4.6.1.  Two -Dimension.!  I  Channel  I- '  ow 


Consider  i  hanne!  >f  width  to,  wi‘ 
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given  by 


and  with  tne  origin  of  foe  y-coordi nates  at  foe  bottom  wa .  i .  me 


asymptotic  analysis  of  the  time-mean  flow  for  the  lower  half  'the 


upper  half  is,  obviously,  symmet' i cal )  provides  results  completi 


analogous  to  those  for  a  boundary  layer,  presented  in  section  4.3,  if 


the  length  scale  A  is  replaced  by  h.  This  time,  however,  the  length 


scale  in  the  longitudinal  1i recti  on  is  infinitely  large  and  need  not 


be  considered.  In  turn,  it  is  fairly  easy  to  show  that  the 


equivalent  to  Eq.  (4.15)  is 
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(4.103) 


The  second  tern  on  the  right  hand  side  of  Eq.  (4.103)  is  part  of  the 


pressure-gradient  contribution.  Again,  the  discussion  at  the 


beginning  of  section  4,4  applies,  with  u  now  defined  as 


(4.104) 


The  outer  scaling  is,  in  analogy  with  Eq.  (4.1M) . 
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Introducing  F.q.  (4.105)  into  l).  (4.103)  and  taking  t  he  limit 


R*  =  u*h/v  ►  =°,  one  gets 


(4. 106) 


which,  obviously,  qi  /“s  Ire]  /  so ! '  ■oreeorv  i  ng  solutions  in  the  form 
of  3  "  v  e  1  oc  1 1  y -d >  .• 1 1-  _  t  1  jv( ' ,  'idii’i;' 
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Equation  (4. 10/)  is  cornu 1  e!  •  /  ana 1  agues  to  Eq.  (4.3/).  beyond  the 
low  frequency  regime  (i.  . ,  for  -t  >>  1),  Lq.  (4.106)  shows 
that  ij  _  0  at  the  lowest  order,  i  ,e. ,  the  outer  layer  or  "core 
region"  undergoes  sleg-fln;  o-,  -i  !  latinos,  as  in  boundary- 1  ayer  flow. 

The  inner  scaling,  on  the  other  hand,  is  the  same  as  in  the 
boundary -layer  case,  i .  n .  ,  it  is  given  by  Eq.  (4. 38) .  To  the  zero 
order,  Eq.  (4.10  3)  tine  :u  ,  -s  t  > 


Eq.  (4.108)  is  •  h  •  ru  mt-TjMi  *  *f  r  (4.41).  ’ho  sol  it  ions  to  this 
equation  will  on  in  ;  in-  r.»r  of  i  "law-  if -fna-we  1  1" ,  ii-'e  that  at  Eq . 
(4.43).  !he  c  .  ..  ,  ...  -)f  s  ( 4.108)  , 

compared  to  '  : .  ;  :  :  •<  •  *  l  n  >  ‘  ?  width  of 


flow  on  the  inner  or  sur.ace"  .  aye  r .  lo“  :  tna ;  jive;  ••  nr,.- 
contribution  only  in  the  low  t  •••..*.;  i.  :.y  r'ayime.  In  fact,  for  :■  > 
the  flow  in  the  core  region  is  slug-like  and,  then,  d  ■  J x/  Jy  •}.  In 

the  high  and  very  high  f requenoy  regimes,  Ji  neeos  to  he  rera'c! 
according  to  Eq.  ( 4.56),  as  in  boundary -layer  flow,  giving  the  sa  le 
solutions,  namely  Eqs.  (4.69)  am  (4.60) ,  respecti  vely. 

The  results  for  the  unsteady  1 ayer,  in  the  interned! ate  a»J  n ' j 1 
frequency  regimes,  are  completely  analogous  to  that  for  the  boundary 
layer  case. 


4.6.2.  Pipe  Flow 

Very  similar  results  are  obtained  for  pipe  flow.  If  the  radius 
of  the  pipe  is  R,  the  centerline  velocity  is  given  by  Eg.  (4.103)  and 
the  origin  of  the  y-axis  is  chosen  at  the  inner  surface  (a  rather 
unusual  choice  for  an  axi symmetric  flow),  then  the  time-mean  flow 
satisfies  the  same  conditions  as  in  boundary-1  ay or  flow  (section 
4.3),  if  A  is  replaced  by  R.  The  equation  equivalent  to  Eq.  (4.103) 
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If  the  outer  scaling  of  E  ;.  (4.106)  is  use  !,  >.  i  th  h  rep :  mm 


;  jy  R . 


the  lowest  order  approximation  of  lq.  (l.ion)  i  u  the  curia  region 
becomes 


lu)  Uj  =jfrn7dn  *  {Tl~n)"  <Tn  [(1~p'  Tl]Vl  {4'110) 

Eq.  (4.110)  is  the  counterpart  of  Eq.  (4.105),  and  gives  also 
solutions  of  the  form  of  Eq.  (4.10/),  with  h  replaced  by  R. 

Likewise,  the  equivalent  to  fq.  (4.10b)  is 
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(4.111) 


Equation  (4.111)  gives  solutions  in  the  form  of  Eq.  (4.43).  Note 
that,  as  in  channel  flow,  the  i  ern  on  the  right  hand  side  of  Eg. 
(4.111)  gives  a  nonnegl igible  contribution  only  in  the  low  frequency 
regime.  In  fact,  for  >>  1  born  dlJ]/dy  and  vani sh  in  the  core 
region.  That  must  he  zero  can  he  easily  proved  by  considering 


sen  es 


expansions  in  powers  of  for  both  0^  and  V j ,  which, 


inserted  into  the  outer  equation,  give  a  set  of  homogeneous  equations 
with  homogeneous  boundary  conditions. 

For  the  unsteady  layer,  finally,  the  equivalent  to  Eq.  (4.65)  is 
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^/here  5  D  ,M  and  the  standard  srali^-'j  ‘-ms  bnon  -r 
provides  solutions  of  the  Lin  ■>*  .  |.  ( 4.  . 
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d.  Eq.  (4.112) 
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VIEW  :i  THE  THEORY 


5.1.  Gener  aj_ 

The  theory  developed  ’  n  the  previ  ojs  chapter  is  now  applied  to 
the  experimental  results.  This  not  only  constitutes  a  check,  hut 
also  allows  the  identification,  within  experimental  error,  of  so  no  of 
the  similarity  laws  to  which  the  theory  leads.  In  addition,  a  clear 
picture  starts  to  emerge  regarding  the  continuous  transition  of  the 
oscillatory  velocity  profiles  from  the  low  to  the  very  high  frequency 
regime.  Some  of  the  concepts  and  ideas  provided  by  the  theory  are 
also  used  to  generalize  its  application  to  other  flow  quantities. 

Table  5.1  shows  the  most  relevant  parameters  for  trie  present 
data  at  each  station  and  for  both  the  oscillating  frequencies.  It  is 
observed  that  the  majority  of  the  situations  correspond  to  the 
intermediate  frequency  regime.  The  remaining  ones  are  in  the 
transition  range  between  the  low  and  intermediate  frequency 
regimes.  Therefore,  they  are  expected  to  have,  to  some  extent,  the 
main  characteri st i cs  of  both  regimes. 


5.2.  Ve  1  ocj ty_ 

In  this  section,  in  addition  to  lb"  present  data,  some 
experimental  results  f  ran  other  studies  on  i'ouudary -1  ayer,  c.h  uel  . 


and  pipe  flows  are  examined.  fni:  w i  .  1  not  only  demonstrate  the 
generality  of  the  theory,  but  will  also  provide  a  more  comprehensive 
picture  of  the  effect  of  unsteadiness  on  wall-bounded  turbulent  shear 
flows.  Table  5.2  presents  a  list,  of  these  additional  experiments, 
together  with  their  most  '*elev.jut.  parameters. 

The  cases  in  the  low  or  trails i tiona  1  low-intermediate  frequency 
regime  must  be  expected  to  exhibit  the  "‘logarithmic  !  aw" ,  Eg.  (4.53) 
or  Eqs .  (4.74)  and  (4.75).  figure  5.1  presents  the  in-phase 

oscillatory  component,  1J^,  in  the  inner  t.oordi  nates ,  for  experiments 
Ml,  M5,  CA,  T1  and  01.  Data  correspu-'di  ng  to  intermediate  and 
intermediate/high  frequency  regimes  (exner iments  02,  Ji  and  04)  are 
also  shown.  A  line  corresponding  to  the  universal  logarithmic  law 
for  the  mean  velocity  profile  in  steady  flow  is  drawn.  Its  slope  is 
used  as  a  reference.  If  is  read i !y  observed  that,  in  all  the  low  and 
low-intermediate  frequency  regime  experiments,  there  is  a  range  in 
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constant  value  behaves  li-e  a  saMI-  -joint.  These  rest:  ;  >  snow  tr. , 
like  the  time-mean  flow,  the  «v;oi  |  i a  lory  motion  in  Lite  -trier  ’..y-r 
depends  strongly  on  the  part icul ar  c' ow  conditions. 

It  is  possible  to  obtain  Some  lore  i  ■->  f  o r' ;  o  c.- C  i  .in  fro  t  t 
foregoing  figures.  Comparing  ■--* <p  >r inent s  Ml,  MS,  CA  an  i  1  i  i  t  •  it. 

5.1,  it  can  be  concluded  that  the  additive  c  instant  io  the 

logarithmic  law  [Bu  in  Eq,  (4.7b)]  increases  wirii  ....  The  same 

effect  is  produced  by  the  . .me -mean  adverse  pressure  gradient,  as 

seen  from  the  data  for  experiment  11.  It  is  more  difficult  to 

identify  trends  in  the  case  of  the  constant  associated  with  the  out- 

of-phase  component.  This  is  because,  since  it  is  small  in  the  zero 

or  small  time-mean  pressure  gradient  cases,  it  is  subjected  to  larger 
experimental  errors.  However,  since  it  has  to  vanish  for  ^>0 

(quasi  -steady  flow)  and  for  (slug  flow,  except  in  a  very  thin 
Stokes  layer)  a  nonmonotonic  trend  has  to  be  expected,  as  is,  in 
fact,  suggested  by  the  data.  Note  that  the  adverse  pressure 

gradient,  in  turn,  causes  the  constant  to  become  positive  (which 
men;  ‘hit  there  is  a  phase  lag). 

.« I  r.;.i  ly  petitioned,  figures  h.  I  and  5,  v  also  present  the 

y.e  experiments  .12,  J3,  and  J4.  The/  art.  particularly  us  el  u  I 
,  -  el.it  i  /<>ly  accurate  data  extending  to  y  1  were  obtained  in 

i ■  ion*  s .  It  i  *;  observed  ttia'  expen  i  .s  :  ind  Ah  show 
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experiments  J2  and  J4  show  that,  with  increasing  u,  a  slug-flow 
region  develops  in  the  outer  layer  beginning  from  the  outer  edge. 
The  out-of-phase  component,  in  turn,  is  reduced  in  the  inner  layer, 
relative  to  the  lower  frequency  cases  J1  and  J3.  In  fact,  higher 
frequencies  cause  the  local  velocity  to  lead  the  free-stream 

velocity.  This  trend  is  opposite  to  that  produced  by  the  time-mean 
pressure  gradient,  which  tends  to  produce  a  phase  lag  in  the  local 

velocity.  It  is  also  important  to  note  that,  at  a  given  time-mean 
adverse  pressure  gradient  (experiments  J^,  J2  or  J3,  J4),  there  is  no 
significant  change  in  the  distribution  of  in  the  inner  layer. 

This  supports  the  theory  of  Chapter  4  that  the  behavior  of  in  the 
inner  layer  is  quasi-steady  in  the  low  to  intermediate  frequency 
range. 

The  previous  discussion  permits  one  to  construct  a  general, 
though  still  speculative,  picture  of  how  the  oscillatory  components 

behave  when  the  frequency  parameter  w  increases.  Figure  5.3  refers 

to  the  case  of  a  zero  time-mean  pressure  gradient  boundary  layer.  It 

is  observed  that  there  is  a  continuous  transition  between  the  "log  + 

wake"  profile  for  the  low-frequency  regime  and  the  Stokes  (viscous) 
solution  for  the  very  high-frequency  regime. 

According  to  the  theory,  in  the  i ntermedi ate-f requency  regime 
there  must  exist  an  outer  region  of  slug  flow.  That  this  is  in  fact 
the  case  for  the  present  experimental  results  has  already  been  shown 
in  Chapter  3  (see  Figs.  3.21  and  3.22).  The  same  is  true  for  the 
other  experiments  discussed  in  this  section,  as  already  mentioned 
with  reference  to  experiments  J2  and  J4.  Immediately  next  to  this 


layer,  and  closer  to  the  wall,  an  intermediate  unsteady  layer 
develops  where  the  "velocity  defect  law",  Eq.  (4.66)  or  Eqs.  (4.82) 


and  (4.83),  must  be  satisfied.  Figures  5.4  and  5.5  show  the  in-phase 
and  out-of-phase  velocity  components,  respectively,  for  all  the 
present  experimental  data,  plotted  in  the  coordinates  correspondi ng 
to  the  unsteady  layer  [see  Eq.  (4.62)].  The  error  in  the  asymptotic 
results  is  of  the  order  (In  R*)"*.  Admitting  an  error  of  ±0.5  (see 
Table  5.1),  we  conclude  that  the  collapse  is  reasonably  good  for  the 
truly  intermediate  frequency  cases,  in  agreement  with  the  theory. 
Deviations  occur  for  the  transitional  low/i ntermediate  frequency 
cases,  specially  for  the  out-of-phase  component  and  closer  to  the 
outer  edge  of  the  unsteady  layer.  In  these  cases,  the  slug-flow 
region  has  nearly  disappeared,  so  the  unsteady  layer  extends 
practically  up  to  the  edge  of  the  boundary  layer.  Hence,  the 
deviations  must  result  from  the  interaction  with  the  free  stream.  In 
all  the  cases,  the  scatter  for  the  out-of-phase  component  is  much 
larger  than  for  the  in-phase  component.  This  is,  as  already 
mentioned,  due  to  the  larger  experimental  errors  in  the  measurement 
of  Uj2*  It  is  significant  to  note  that  the  location  of  the  peak  in 
the  in-phase  component,  which  is  a  measure  of  the  extent  of  the 
unsteady  layer,  corresponds  to  y  =  1. 

Figures  5.6  and  5.7  present  the  results  for  experiments  T2,  J2, 
J4,  AR  and  BK  plotted  in  the  "unsteady-layer  coordinates".  A  typical 
set  of  data  from  the  present  experiments  is  also  shown  in  the  figure 
for  comparison.  It  can  be  seen,  again,  that  this  normalization 
brings  the  data  closer  to  one  another,  even  though  considerable 
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deviations  still  remain.  However,  these  should  be  interpreted  with 
caution  due  to  the  following  reasons: 

(i )  Except  for  the  data  set  of  experiment  T2  and,  of  course, 
the  present  ones,  all  of  the  data  were  obtained  from 
published  graphs. 

(ii)  The  data  of  experiments  J2,  AR  and  BK,  though  closer  to 
satisfy  the  smal 1 -amp! itude  assumption,  may  be  subjected 
to  larger  experimental  errors  because,  precisely,  of  the 
very  low  amplitudes  of  oscillation  (see  Table  5.2). 

(iii)  In  using  the  data  set  J2  and  J4,  the  local  boundary- 
layer  thickness  used  by  the  originators  to  normalize  the 
y-coordinates  in  their  plots  was  not  available  and  was 
therefore  estimated. 

There  is  a  relatively  good  collapse  of  the  results  for  the  in-phase 
component  in  the  near-wall  region  (see  Fig.  5.5).  In  the  outer  part, 
some  of  the  relative  departures  may  be  genuine  and  linked  to  the 
different  characteristics  of  the  outer  flow.  Similar  comments  can  be 
made  about  the  out-of-phase  component  (Fig.  5.7)  which,  again,  show 
significant  scatter  in  the  near-wall  region. 


5.3.  Other  Flow  Properties 

Though  the  theory  has  been  developed  only  for  toe  oscillatory 
velocity  components,  the  ideas  can  be  extended  to  analyze  other  flow 
properties  also.  If  the  in-phase  and  out-of -phase  components  of  the 

9  9 

turbulence  properties,  namely  <u  >,  <v  >,  and  <uv>,  are  defined  in 

complete  analogy  with  those  for  the  velocity  [see  Eq.  (4.71)],  they 
2 

must  scale  with  u*.  To  ue  consistent  with  eariler  definitions  (see 

2 

chapter  3),  they  are  normalized  with  2u*.  Their  distributions  in  the 
unsteady -layer  coordinates,  for  the  present  experiments,  are  shown  in 
Figs.  5.8  to  5.13.  With  the  same  admitted  error  as  for  the  velocity 
(±0.5),  the  data  for  the  in-phase  component  collapse  fairly  well  for 

p 

all  the  quantities,  except  for  <u  >  in  the  outer  part  of  the  unsteady 
layer.  Deviations  occur  there  for  the  transitional  1 ow/i ntermediate 
frequency  cases  (see  Fig.  5.8).  As  in  the  case  of  the  velocity,  this 
must  be  attributed  to  interactions  with  the  free  stream,  including 
some  feeding  of  free-stream  turbulence  into  the  boundary  layer.  Less 
satisfactory  is  the  performance  of  the  out-of -phase  components  (Figs. 
5.11  to  5.13),  especially  in  the  near-wall  region.  This  can  be  due 
to  a  higher  sensitivity  of  the  turbulent  quantities  to  the  value 
of  a.  Note  that  the  best  data  collapse  is  for  <uv>.  The  agreement 

9 

improves  when  moving  outwards,  except  for  <u  >,  for  which  the  same 
considerations  as  for  the  in-phase  component  apply. 

As  in  the  case  of  the  turbulence  properties,  the  in-phase  and 
out-of-phase  components  of  the  wall  shear  stress  <tw>  (obtained  from 
<Cf>)  must  scale  with  ^2  Figure  5.14  presents  the  amplitude  of 

oscillation  of  the  wall  shear  stress  as  a  function  of  the  frequency 
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parameter  Included  are  the  measurements  at  all  stations  and  all 
frequencies.  They  merge  satisfactorily  within  the  admitted  error, 
+0.08  (obtained  from  that  for  [<Cf>];  see  Table  2.1).  It  is  observed 
that  the  amplicude  remains  fairly  constant  with  Z.  There  possibly 
exists  a  mild  minimum  around  w  -  30  to  35.  The  variation  of  the 
phase  of  <tw>  with  frequency  is  shown  in  Fig.  5.15.  Again,  there  is 
a  satisfactory  correlation  of  the  results.  An  approximately  linear 
trend  is  observed  in  the  variation  of  the  phase  angle  with  w. 

Finally,  the  in-phase  and  out-of-phase  components  of  the  eddy 

viscosity,  obtained  from  the  measurements  at  station  5,  and  at  both 

the  frequencies,  are  plotted  in  Figs.  5.15  and  5.16.  They  are 
2 

normalized  by  u*/a),  as  suggested  by  dimensional  analysis.  A 
relatively  large  uncertainty  must  be  accepted  because  of  the  manner 
in  which  they  are  calculated.  For  example,  they  are  very  sensitive 

to  the  values  of  al^/ay,  which  are  obtained  from  a  piecewise  fit  to 

vs.  y  (see  chapter  3).  The  collapse  for  the  two  sets  of  data  is 

fairly  good  in  the  near-wall  region,  especially  for  the  in-phase 
component.  Large  differences  are  observed  in  the  outer  part,  which 
are  believed  to  be  associated  with  the  above-mentioned  calculation 
errors.  The  results  corresponding  to  the  lower-f requency  data  set 
are  considered  to  be  more  reliable  in  this  respect.  The  data, 
however,  seem  to  confirn  that  neither  a  quasi-steady  nor  a 

frozen  eddy-viscosity  model  (v  =  v  =  °)  is  realistic  at  either 
of  these  frequencies. 


5.4.  Time-Mean  Flow 


As  discussed  in  chapter  3,  small  but  significant  •. 

effects  are  found  in  the  outer  distribution  of  the  time-mean 
quan  t  i  ties,  specially  the  turbulence  properties.  Though  tile  theory 
developed  does  not  deal  with  the  time-mean  flow  (in  fact,  it  assumes 
that  it  is  equal  to  the  steady  flow  at  the  time-mean  free-stream 
velocity),  the  frequency  parameters  derived  from  the  analysis  can  be 
used  to  correlate  the  observed  departures. 

The  strongest  frequency  effects  were  observed  in  experiments  M4, 
M5,  M6  and  M7.  Hence,  it  can  be  said  that  frequency  effects  on  the 
time-mean  quantities,  in  the  outer  layer,  can  be  expected  in  the 
range  8  <  Z  <  20.  For  w  <  8,  frequency  effects  are  too  weak.  On  the 
other  hand,  for  Z  >  20,  frequency  effects,  if  present,  must  be 
confined  to  the  near-wall  region.  In  the  outer  layer  the  turbulence 
distribution  is  then  frozen,  so  the  time-averaged  flow  properties 
approach  those  in  steady-flow  at  the  time-mean  velocity. 
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Figure  5.1.  In-phase  velocity  component  for  the  low  and 
low/intermediate  frequency  regime.  A,  experi¬ 
ment  Ml:  □  .  M5:  *.  CA;  0  .  T1 ;  y,  Jl;  '  ,  J2; 

+  ,  J3;  x,  J4; - ,  universal  logarithmic  law. 
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Figure  5.2.  Out-of-phase  velocity  component  for  the  low 
and  low/intermediate  frequency  regime. 
Symbols  as  in  Fig.  5.1 . 
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Figure  5.9.  In-phase  component  for  <v  >  in  unsteady  layer 
coordinates.  Symbols  as  in  Fig.  5.4. 


Figure  5.16.  In-phase  component  of  <v.>  in  unsteady  layer 
coordinates.  □  ,  experiment  M5;  0,  M10. 


CHAPTER  6 


CONCLUSIONS 


The  following  are  the  main  findings  of  the  present 

investigation: 

(i)  Spatial  history  effects  on  the  time-mean  properties  have 
been  identified  and  interpreted  for  a  periodic,  turbulent  boundary 
layer  developing  in  zero  time-mean  pressure  gradient.  It  has  been 
found  that  the  oscillation  accelerates  the  development  of  the  outer 
flow,  while  it  has  little  effect  on  the  inner  flow.  The  modification 
of  the  outer  flow  can  be  interpreted  as  a  shift  in  the  virtual  origin 
of  the  turbulent  boundary  layer.  For  large  Reynolds  numbers,  the 
time-mean  properties  tend  to  approach  their  steady  distributions  at 
the  time-mean  free-stream  velocity.  The  small  residual  differences 
between  the  two  can  be  considered  to  be  'frequency  effects'. 

(ii)  Frequency  effects  on  the  time-mean  properties  have  also 
been  identified  and  correlated.  They  produce  a  small  positive 
departure  in  the  outer  layer,  relative  to  the  steady-flow 
distribution  at  the  same  Reynolds  number,  for  8  <  Z  <  20. 

(iii)  The  behavior  of  the  oscillatory  components  of  the  motion 
has  been  seen  to  confirm  what  was  previously  known.  In  particular, 
many  of  the  conclusions  of  Ramaprian  and  Tu  (1982)  on  periodic  pipe 
flows  apply  equally  well  to  unsteady  boundary  layers.  The  thickness 
of  the  unsteady  layer,  i.e.  the  layer  over  which  unsteady  viscous 


effects  are  significant,  has  been  characterized  by  the 
quantity  u */w,  in  agreement  with  Ramaprian  and  Tu  (1982).  Phase 
differences  of  about  360  degrees  in  the  oscillatory  components  of  the 
turbulence  properties,  across  the  unsteady  layer,  have  been 
measured.  Eddy-vi scosity  calculations,  obtained  from  the 

measurements,  have  shown  that  neither  a  quasi -steady  nor  a  'frozen' 
turbulence  model  is  appropriate  for  unsteady  boundary  layers  at  the 
present  frequencies. 

(iv)  Wal 1 -shear-stress  measurements  have  been  made  using  a 
flush-mounted  heat-flux  gage.  The  time-mean  wall  shear  stress  has 
been  found  to  decrease  slightly  with  increasing  frequency,  a  trend 
opposite  to  that  observed  in  pipe  flow  [Ramaprian  and  Tu  (1982)]. 
The  amplitude  of  oscillation  is  close  to  the  quasi-steady  value  at 
all  the  frequencies.  Larger  phase  leads  (of  over  30  degrees)  than 
those  reported  by  Ramaprian  and  Tu  (1982)  have  been  observed  ir.  fhe 
boundary  layer  at  the  higher  frequencies. 

(v)  An  asymptotic  theory,  for  large  Reynolds  numbers,  has  been 
developed  for  the  oscillatory  motion.  It  is  valid  for  both  boundary 
layers  at  zero  and  adverse  time-mean  pressure  gradients  and  fully 
developed  pipe  and  channel  flow,  and  has  been  successfully  applied  to 
the  present  and  previous  available  experimental  information. 

(vi )  The  theory  identified  the  frequency  parameters  u  =  uA/u* 
and  a  =  u*/wv,  as  the  appropriate  parameters  to  characterize  the 
oscillatory  motion  at  large  Reynolds  numbers.  The  first  one  is 
analogous  to  that  introduced  by  Ramaprian  and  Tu  (1932)  based  on  eddy 
vi scosi ty  arguments. 


(viii)  Four  frequency  regimes  have  been  identified  (excluding 
the  quasi -steady  regime): 

*  Low  frequency  regime:  Z  <  1  (a  >  R*) 

~  1/9  1/9 

*  Intermediate  frequency  regime:  1  <<  u  <  R*  (R*  <  a  <<  R^; 

*  High  frequency  regime:  R <<  Z  <<  R*  (1  <<  a  <<  r]/2) 

*  Very  high  frequency  regime:  Z  >  R*  (a  <  1) 

Similarity  laws  have  been  identified  for  each  one  of  these  frequency 
regimes.  These  laws  are  the  extension  to  unsteady  flows  of  the  well 
known  "1 aw-of -the-wal  1 ",  "velocity-defect  law"  and  "logarithmic  law" 
for  steady  flows,  and  are  summarized  in  section  4.3.5. 
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